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ABSTRACT 

This  study  develops  simulation  models  for  the  evaluation  of  the  traffic  performance 
at  a  three  way  intersection  with  no  signal.  In  particular,  the  models  were  designed  to  test 
and  to  evaluate  the  wartime  road  capacity  of  the  current  and  a  proposed  road  system  in 
the  Korean  theater.  The  models  describe  how  traffic  at  the  road  intersection  will  vary 
with  changes  to  the  input  variables  and  priority  rules.  An  analysis  is  performed  to  com- 
pare performance  measures  of  the  intersection  given  various  traffic  conditions.  Specif- 
ically, the  analysis  will  enable  an  efficient  modification  of  the  roads  or  an  appropriate 
road  usage  plan  in  wartime. 
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I.      INTRODUCTION 

A.      BACKGROUND 

The  wartime  traffic  on  the  roads  in  the  southern  theater  of  the  Republic  of  Korea 
needs  special  attention,  especially  the  traffic  involved  in  a  war  which  is  very  dynamic, 
and  in  particular  a  war  which  involves  a  sudden  attack  from  the  north  and  the  necessity 
of  falling  back  to  the  south.  There  are  several  main  avenues  leading  from  the  front,  from 
which  many  streets  to  the  east  or  west  are  derived,  so  we  have  many  three  way  inter- 
sections with  undeveloped  roads.  The  traffic  on  these  roads  is  not  very  heavy  in  peace 
time,  but  it  appears  as  if  the  limited  road  capacity  may  restrict  the  timely  movements 
of  forces  in  a  time  of  emergency. 

The  study  of  traffic  in  this  specific  area  has  never  been  done  before  because  under 
the  current  strategy  no  forces  are  allowed  to  withdraw,  with  the  exception  of  supply 
forces  and  forces  needing  to  be  replaced.  We  assume  that  in  many  cases  all  the  forces 
will  have  to  move  simultaneously  or  in  a  certain  sequence.  For  instance,  the  enemy  may 
select  battlefields  which  are  quite  different  from  those  we  expect.  Also,  the  deployed 
forces  may  move  forward  or  backward  simultaneously  as  the  situation  changes.  It  is 
stronalv  believed  that  traffic  congestion  will  occur  in  the  assumed  contingencies,  and 
therefore,  it  is  necessary-  to  have  some  idea  of  the  dynamics  involved.  A  particular 
question  which  could  be  asked  is,  under  what  condition  will  the  traffic  moving  onto  the 
main  road  from  the  side  road  be  unstable.  That  is,  when  will  the  queue  of  traffiic  continue 
to  grow  in  an  unstable  manner.  This  thesis  presents  a  ver\'  Umited  approach  to  the 
problem  of  the  dynamics  offeree  withdrawal  from  the  front  or  a  deployment  of  reserve 
forces  to  the  front,  but  it  is  an  important  first  step  to  the  understanding  and  solution 
of  the  whole  problem  of  traffiic  control  and  flow  in  the  Korean  theater  of  operations. 


B.      DESCRIPTION  OF  THE  PROBLEM 

The  purpose  of  building  a  stochastic  simulation  model  is  to  describe  traffiic  per- 
formance at  a  three  way  intersection.  This  type  of  intersection  is  a  principal  capacity 
limiter  and  main  source  of  delay  of  the  current  road  system  in  the  Korean  theater.  At 
first,  an  analytical  model  for  the  current  road  system  under  special  simplifying  assump- 
tions is  developed,  including  a  transition  diagram  and  brief  solution  procedure.   We  will 


assume  that  a  unit,  which  possibly  consists  of  many  vehicles,  is  treated  as  a  single  cus- 
tomer as  in  general  queueing  problems.  Each  unit's  transition  time  is  equated  with  ser- 
vice time  in  the  sense  of  a  classical  queueing  situation,  [Refs.  1,2:  pp.  13-20].  The 
interarrival  and  service  times  of  a  unit  are  assumed  to  have  certain  respective  distrib- 
utions, which  will  be  discussed  in  detail  later.  Further,  the  interaction  of  different  units 
at  the  intersection  will  be  controlled  by  traffic  policies  to  be  described  in  Chapter  II. 

Due  to  difficulties  encountered  in  the  analytical  approach,  two  simulation  models 
are  developed  for  the  current  and  a  proposed  road  system.  We  are  concerned  with 
finding  how  the  traffic  performance  at  the  intersection  varies  with  changes  to  the  arrival 
and  service  time  distributions.  We  consider  performance  in  terms  of  the  average  waiting 
time  of  the  units  and  the  average  number  of  units  waiting  in  each  direction  of  the  inter- 
section. 

The  intersection  system  must  be  thoroughly  understood  because  blockages  at  an 
intersection  may  cause  serious  problems  in  conducting  operations  in  an  area  where  units 
are  deployed  with  high  density  relative  to  the  road  capacity.  The  units  moving  or  wait- 
ing because  of  blockages  are  very  vulnurable  to  enemy  ambushes  since  they  are  not  fully 
ready  for  engagements. 


C.      SCOPE  OF  THE  THESIS 

In  attempting  to  achieve  the  goals  of  this  study,  an  analytical  model  is  developed 
up  to  the  transition  diagram,  together  with  the  general  solution  procedure  outlined  in 
Chapter  II.  It  is  found  that  the  simulation  model  is  needed  to  solve  the  problems  be- 
cause of  the  difficulties  in  handhng  various  inputs  and  because  of  the  limited  solution 
capability  available  for  the  analytical  model. 

Two  simulation  models  are  built  with  SIMSCRIPT  II. 5  in  order  to  conduct  the 
simulation  experiments.  The  first  one  is  based  on  the  current  road  system  and  the  other 
is  for  a  proposed  modification  to  the  system.  The  two  models  are  presented  in  Chapter 
III  and  Chapter  IV.  Figure  1  and  Figure  15  depict  the  current  and  proposed  road  sys- 
tems respectively,  with  possible  movement  routes. 

Chapter  V  will  analyze  and  compare  the  results  of  the  two  simulation  models.  The 
analysis  is  done  by  exporting  output  files  of  SIMSCRIPT  runs  to  the  GRAFSTAT 
program  on  the  IBM  3033AP  main  frame.  This  allows  one  to  graphically  examine  the 
models  for  stability  of  the  performance  measurements  and  to  examine  the  results  in  the 


stable  cases  as  a  function  of  several  input  factors.    Finally,  conclusions  and  recommen- 
dations are  presented  in  Chapter  VI. 


II.      PROBLEM  STATEMENT 

As  mentioned  before,  we  can  apply  this  model  to  the  wartime  traffic  in  a  theater  of 
operations  where  the  road  capacity  is  limited  in  comparison  with  the  potential  road  us- 
ers. We  assume  that  there  are  no  traffic  signals  and  no  military  police  for  traffic  control 
and  assume  services  are  controlled  by  service  policies  we  set  up  for  the  simulation  study. 
It  is  the  fundamental  road  capacity  and  traffic  performance  at  the  intersection  that  we 
want  to  examine,  therefore,  we  must  fully  understand  the  following  assumptions  to 
properly  grasp  the  specific  features  and  logic  of  the  model.  A  current  intersection  is 
presented  in  Figure  1. 

A.      ASSUMPTIONS 

1.  Units 

A  unit  consists  of  a  group  of  vehicles,  the  number  of  which  depends  on  the  or- 
ganization of  forces  and  the  tactics  of  a  war.  All  the  vehicles  in  the  unit  will  be  treated 
as  one  customer  at  the  intersection  (general  queueing  facility)  because  the  movement  of 
a  combat  unit  such  as  a  company  or  a  battery  must  be  conducted  under  tight  command 
and  control  to  minimize  confusion  and  loss  of  vehicles.  Thus  we  assume  that  no  unit 
or  vehicle  can  intervene  between  another  unit's  vehicles. 

2.  Arrivals  and  Queues 

Poisson  arrivals  are  assumed  for  the  analytical  model.  Each  arrival  from  the 
three  directions  is  filed  in  an  appropriate  queue  if  it  is  not  its  turn  for  service.  Our 
intersection  queueing  facility  has  three  queues;  a  north,  south  and  east  queue  and  there 
is  no  interaction  between  vehicles  in  each  of  the  queues,  that  is,  a  First  Come  First 
Served  principle  is  effective  within  each  direction.  Actually  the  time  of  arrival  is  that  of 
the  unit  to  the  end  of  the  hne,  potentially  some  distance  away  from  the  intersection. 
Since  units  move  up  to  the  intersection  in  close  order  as  units  which  arrived  previously 
are  served,  the  distance  and  stop  and  start  effects  of  the  queue  at  the  intersection  are 
ignored.  In  fact,  the  waiting  time  of  a  unit  is  proportional  to  its  distance  from  the 
intersection.  This  could  be  exploited  in  the  analysis  but  will  be  ignored. 

3.  Route  Selection 

Possible  routes  for  a  unit  from  each  direction  are  shown  by  arrows  in  Figure  I . 
One  of  two  routes  will  be  selected  for  a  unit  from  each  direction  and  this  selection  is 
made  when  a  unit  becomes  the  first  unit  in  the  queue  and  no  unit  is  presently  served 
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route2 
East  route  1 


North  route2 
South  route2 


South  route! 


Figure  1.  A  current  Intersection.:  Traffic  coming  from  the  north  either  contin- 
ues south  or  turns  east  if  it  is  not  blocked.  Similarly  traffic  from  the 
south  can  continue  north  or  turn  east  if  it  is  not  blocked  and  traffic  from 
the  east  turns  either  north  or  south.  Priorities  for  service  on  the  routes 
will  be  discussed  in  detail  in  Chapter  II. 

from  the  direction  or  when  an  arriving  unit  sees  no  unit  in  the  queue.  The  probability 
used  to  decide  a  route  is  an  input  variable  in  the  model  and  a  uniform  (O.I)  random 
variable  U  will  be  drawn  to  determine  the  route  selection.  If  U  is  less  than  the  proba- 
bihty  P  of  taking  route  1,  route  1  is  selected,  otherwise  route2  is  selected.  The  probabihties 


of  taking  north  route  1,  south  route  1,  and  east  route  1  may  be  difTerent  from  each  other 
and  are  user  input  parameters. 

4.  Service 

Service  begins  when  the  first  vehicle  of  a  unit  starts  from  the  intersection  and 
ends  upon  the  passage  of  the  last  vehicle.  Once  served  a  unit  is  assumed  to  leave  the 
system,  even  though,  in  reality,  it  may  significantly  affect  the  service  conditions  of  follow 
on  units.  Service  times  depend  on  many  factors,  such  as  the  length  of  a  unit,  road  con- 
ditions, weather,  number  of  refugees  and  so  on.  Though  many  of  the  factors  should  be 
studied  separately,  they  are  assumed  to  be  compounded  into  a  single  random  variable 
used  for  service  time  in  this  model.  Our  main  interest  in  this  thesis  is  the  interactions 
between  units  at  the  intersection.  We  assume  the  service  times  for  the  six  routes  are 
independent  and  exponentially  distributed  in  the  analytical  model.  This  is  done  for  an- 
alytical tractability  and  is  probably  not  physically  valid.  More  general  and  probably 
more  realistic  assumptions  are  made  in  the  simulation  models. 

5.  Service  Policies  for  INIodel  1 

In  Model  1,  service  for  a  unit  starts  immediately  upon  arrival  when  there  is  no 
blocking  to  the  selected  route,  otherwise  it  is  filed  in  the  queue.  There  are  two  possible 
routes  from  each  direction  named  route  1  and  route2  as  shown  by  arrows  in  Figure  1. 
Suppose  that  the  set  (north  route,  south  route,  east  route)  indicates  the  state  of  move- 
ments at  the  intersection  where  each  entity  within  that  set  can  have  value  0,  1  or  2  ac- 
cording to  the  routes  being  used.  Here  0  indicates  that  neither  of  the  two  possible  routes 
from  each  direction  are  used;  otherwise  1  and  2  represent  the  route  number  used.  It  is 
clear  that  all  combinations  of  3  ( =  27)  states  are  not  included  in  the  possible  movement 
states  of  the  model  because  of  the  interactions  between  the  routes.  For  instance,  the 
state  (1,1,2)  never  happens  since  south  route  1  and  east  route2  can  not  be  used  simul- 
taneously. 

In  Model  1,  in  fact  there  are  only  14  states  of  movement  which  are  acceptable: 
(0,0,0),  (0,0,1),  (0,0,2),  (0,1,0),  (0,2,0),  (0,2,1),  (0,2,2),  (1,0,0),  (1,0,2),  (1,1,0),  (1,2,0) 
(1,2,2),  (2,0,0),  and  (2,0,2).  Even  though  these  sets  are  not  themselves  enough  to  study 
the  quantities  of  interest,  the  service  polices  are  based  on  these  possible  movement 
combinations. 


A  movement  on  a  certain  route  may  block  more  than  one  route  to  which  units 
have  to  move,  therefore  the  priorities  for  movement  must  be  specified  according  to  the 
route  for  which  a  service  is  just  completed.  Note  that  we  have  assumed  there  is  no  traffic 
signals  and  no  military  pohce  for  traffic  control.  A  unit  desiring  to  move  in  the  system 
must  wait  until  its  route  is  clear  and  it  has  priority  (in  the  case  of  ties).  We  break  ties 
in  cases  in  which  more  than  one  unit  has  conflicting  priorities  simultaneously.  The 
routes  that  block  movement  and  the  next  priority  for  a  service  given  a  blocking  unit 
which  has  completed  service  are  shown  in  Table  1. 

Table   1.     PRIORITIES  AT  THE  END  OF  SERVICE  OF  A  BLOCKING  ROUTE 
(MODEL  1) 


Route  which  is  about 
to  finish  service 

Next  Priority,  in  order,  at  the  end  of  blocking 

north  route! 
north  route2 
south  route  1 
south  route2 
east  route  1 
east  route2 

east  routel,  (north  routel,2) 
east  routel,  (south  routel, 2),  (north  routel, 2) 
(east  routel, 2),  north  route2,  (south  routel, 2) 

north  route2,  (south  routel, 2) 
(north  routel, 2),  south  routel,  (east  routel, 2) 

south  routel,  (east  routel, 2) 

The  left  columi^  of  each  row  in  Table  1  indicates  the  routes  on  which  service 
was  completed  and  the  right  column  of  each  row  identifies  all  possible  routes  which  can 
be  selected  by  units  waiting  to  pass  when  service  was  completed  by  the  route  denoted  in 
the  left  column.  Priorities  for  next  service  are  given  to  the  right  column  routes  in  the 
order  Usted.  We  set  up  the  specific  priorities  for  this  model  by  breaking  ties  and  placing 
the  routes  for  which  service  was  just  completed  last  in  order  of  priority.  Next  priorities 
at  the  end  of  service  on  each  route  will  be  discussed  in  detail  below.  The  routes  in  the 
parenthesis  denote  the  two  possible  candidates  from  the  direction  blocked,  so  only  one 
of  the  two  routes  can  be  selected  for  the  first  unit  in  the  queue.  Table  1  also  shows  that 
a  newly  arriving  unit  desiring  to  move  on  a  route  of  a  left  column  and  row  can  only 


move  if  all  the  right  column  routes  of  the  same  row  are  not  used,  othen\'ise  it  has  to  wait 
for  service  until  next  priority  is  given  to  it. 

After  service  for  a  unit  on  a  route  in  the  left  column  of  Table  1  is  completed,  the 
priority  for  next  service  is  conditioned  on  other  movements  at  the  intersection.  As  an 
example,  the  third  row  could  represent  units  that  have  to  move  to  (east  route!  or  east 
route2),  north  route2,  and  (south  route]  or  south  route2)  and  are  blocked  by  a  unit 
presently  using  the  left  column  south  route! .  Upon  service  completion  on  south  route! 
the  first  check  is  to  determine  if  there  is  blocking  of  east  route!  or  east  route2.  Let's 
assume  that  a  unit  from  the  east  has  been  waiting  for  movement  to  the  north  (east 
route2  in  the  first  position  of  the  right  column),  the  unit  can  move  immediately  because 
east  route2  is  always  clear  right  after  the  service  completion  of  south  route!.  Next, 
service  for  north  route2  (second  position  in  order)  will  be  considered  and  a  unit  waiting 
to  go  from  the  north  to  the  east  (north  route2)  can  also  move  immediately  because  there 
is  no  blocking  to  that  route  at  the  end  of  south  route  !'s  service  regardless  of  the  service 
on  east  route2.  Note  that  east  route2,  which  has  the  first  priority  for  service,  does  not 
block  north  route2.  Finally,  service  on  the  south  routes  (third  position  in  priority  order) 
will  be  considered  and  movement  to  either  south  route!  or  south  route2  will  be  condi- 
tioned on  the  status  of  east  route2  and  north  route2.  If  we  had  assumed  east  route!  (the 
other  route  in  the  first  parenthesis  of  the  right  column  in  the  third  row)  was  waiting  to 
move  instead  of  east  route2,  we  must  check  to  determine  if  there  is  a  unit  on  north 
route!  upon  the  service  completion  of  south  route!.  If  north  route!  is  active  this  route 
also  blocks  east  route!.  If  there  is  no  unit  on  north  route!,  east  route!  can  be  used, 
otherwise  service  on  the  next  route  in  the  order  of  priority,  north  route2,  must  be  con- 
sidered. In  other  words,  even  though  a  unit  has  priority  upon  a  route's  service  com- 
pletion, the  unit  loses  priority  if  it  is  blocked  by  route  presently  in  service. 

In  Table  2  through  Table  7  we  tabulate  cases  of  the  priority  problems.  The 
triple  in  the  row  designator  of  each  Table  describes  the  routes  selected  by  units  which 
are  specifically  blocked  by  the  route  that  is  about  to  finish  service.  The  first  element 
refers  to  the  traffic  from  the  north,  the  second  refers  to  traffic  from  the  south  and  the 
third  refers  to  traffic  from  the  east.  In  each  case  the  element  can  have  value  0,  1  or  2, 
where  0  indicates  no  traffic  waiting  from  that  direction,  1  indicates  traffic  waiting  for 
route!  and  2  indicates  that  trafi'ic  is  waiting  for  route2.  The  triples  in  the  column  des- 
ignator represent  all  possible  movement  states  when  the  route  of  interest,  which  is 
specified  in  the  Table  title,  is  about  to  finish  service.    P'inally,  NI,  N2,  SI,  S2,  EI,  and 


E2  represent  north  route  1,2,  south  route  1,2,  and  east  route  1,2  respectively.  A  route 
identified  at  the  intersection  of  a  row  and  column  has  priority  upon  the  service  com- 
pletion of  the  specified  unit  in  the  title  of  the  Table.  N/A  indicates  that  the  combination 
of  the  row  and  column  condition  is  not  feasible  in  this  Table.  A  0  at  the  intersection  of 
a  row  and  column  indicates  that  no  additional  route  can  be  used  immediately  after  ser- 
vice completion  of  the  unit  specified  in  the  Table  title. 


Table  2.     SERVICE  PRIORITIES  AT  THE  END  OF  A  NORTH  UNIT  MOVING 
SOUTH  (denoted  by  north  route  1) 


N    S    E 

1     0    0 

1     0    2 

1     1     0 

1     2    0 

1     2    2 

0  0     I 

1  0    0 

1  0     1 

2  0    0 
2    0     1 

El 
Nl 
El 
N2 
El 

N/A 
Nl 

N/A 
N2 

N  A 

0 

NT 

Nl 

0 

0 

El 
Nl 
El 
0 
El 

N/A 
NT 

N/A 
0 

N/A 

In  the  last  row  of  Table  2  the  entry  (2  0  1)  in  the  left  column  means  that  when  service 
on  north  route  1  is  completed,  north  route2  is  selected  for  the  next  north  unit  (the  first 
2),  there  is  no  one  waiting  from  the  south  direction  (the  second  0)  and  someone  is 
waiting  to  go  from  the  east  to  the  south  (the  last  1 ).    The  priorities  for  movement  of 
waiting  units  must  be  conditioned  on  the  present  situation  at  the  intersection,  and 
this  is  the  column  indicator.  Thus  while  a  unit  from  the  north  is  moving  south,  a  unit 
from  the  south  can  be  moving  north  (south  route),  designated  by  the  second  1  in  the 
header  of  column  3)  if  no  unit  from  the  east  is  turning  north  (east  routc2).  and  can 
move  east  at  any  time  (south  route2).  Therefore  there  are  only  five  movement  states 
possible  when  north  route!  is  about  to  finish  service.  The  priorities  for  next  service 
given  the  row  entry  are  dependent  on  which  one  of  the  five  movement  states  describes 
the  intersection  just  prior  to  the  completion  of  north  route  1.  Not  all  row  and  column 
combinations  are  feasible  and  are  designated  by  N/A.  Thus  in  row  (2  0  1)  and  column 
(1  0  2),  east  route2  is  being  used  so  no  route  is  selected  for  an  east  unit  (the  1  of  (2 
0  1)  can't  be  selected).  A  0  indicates  that  no  additional  route  can  be  used  after  service 
completion  of  north  route  1.    Thus  for  row  (2  0  0)  and  column  (1  2  2),  the  unit  from 
the  north  can  not  move  east  because  a  unit  from  the  south  is  moving  east  (2  in  the 
second  position  of  the  column  heading). 
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Table  3.  SERVICE  PRIORITIES  AT  THE  END  OF  A  NORTH  UNIT  TURNING 
EAST  (denoted  by  north  route2):  Explanation  of  Table  arrangement  is 
provided  on  pages  8  and  9. 


N    S 

E 

2    0    0 

2    0    2 

While  a  unit  from  the  north  is  turning 
east  (north  route2),  only  east  route2  is 
not  blocked  by  the  unit.  Thus  we  have 

0    0 

1 

El 

N;A 

0     1 

0 

SI 

0^ 

only  two  possible  movement  combina- 

0    1 

1 

El 

N/A 

tions  at  the  intersection  just  prior  to 

0    2 

0 

S2 

S2 

north  route2  finishing  service.    All  com- 

0 2 

1  0 
1     0 

1 

0 

1 

S2  and  El 
M 
El 

N/A 
Nl 

N/A 

binations  of  the  6  routes  except  east 
route2  are  in  the  row  designators,  which 

indicate  the  possible  routes  selected  by 
units  waiting  to  pass  when  service  on 

1     1 

0 

Nl  and  SI 

.vr 

north  route2  is  completed.     A  unit  on 

1     1 

1 

El 

N/A 

east  route2  (the  last  2  in  column  (2  0  2)) 

1     2 

0 

Nl  and  S2 

Nl  and  S2 

still  blocks  a  unit  waiting  for  service  on 

1     2 

1 

El  and  S2 

N/A 

south  routel  when  north  route2's  service 

2    0 
2    0 
2     1 

0 

1 

0 

N2 
El 
SI 

N2 
N/A 
A'2^^ 

is  completed,  thus  south  routel  loses 
priority,  which  is  designated  by  V  .  The 
case  of  W  in  row  (2  1  0)  and  column  (2 
0  2)  indicates  that  the  priority  initially 

2     1 

1 

El 

N/A 

given  to  south  routel  is  transferred  to 

2    2 

0 

S2 

S2 

north  route2  because  of  the  unit  on  east 

2    2 

1 

El  and  S2 

N'A 

route2. 
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Table  4.  SERVICE  PRIORITIES  AT  THE  END  OF  A  SOUTH  UNIT  MOVING 
NORTH  (denoted  by  south  routel):  Explanation  of  Table  arrangement 
is  provided  on  pages  8  and  9. 


N    S 

E 

0     1     0 

1     1     0 

South  routel    blocks  all  routes   except 

0    0 

1 

El 

o' 

0    0 

2 

E2 

E2 

north  routel.  Thus  we  have  two  possible 
movement    combinations   just   prior   to 

0     1 

0 

SI 

SI 

south  routel    finishing   service   and   all 

0     1 

1 

El 

sr 

combinations  of  routes  which  have  been 

0     1 

2 

E2 

E2 

blocked  by  a  unit  on  south  routel  in  the 

0    2 

0 

S2 

S2 

row  designators  to  identify  priorities  at 

0    2 

1 

S2andEl 

52^ 

the  end  of  service  on  south  routel.     A 

0    2 
2    0 
2    0 
2    0 

2 
0 

1 

2 

S2  and  E2 
N2 
El 

N2  and  E2 

S2  and  E2 
N  A 
N  A 
N  A 

unit  from  the  north  moving  south  (the 
first  1  of  column  (1  1  0))  blocks  east 
routel  again  after  service  completion  on 
south  routel.  This  is  desisnated  bv  V  in 
this  Table.    As  an  example,  in  row  (0  2 

2     1 

0 

N2 

N  A 

1)  and  column  (1   1  0),  priority  initially 

2     1 

1 

El 

N/A 

given  to  cast  routel  must  be  transferred 

2     1 

2 

N2  and  E2 

N  A 

to  south  route2  because  north  routel  is 

2    2 

0 

\2 

X  A 

being  used  and  blocks  east  routel  when 

2    2 

1 

El  and  S2 

N  A 

service  on  south  routel  is  completed. 

2    2 

2 

N2  and  E2 

N  A 
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Table  5.     SERVICE  PRIORITIES  AT  THE  END  OF  A  SOUTH  UNIT  TURNING 
EAST  (denoted  by  south  route2) 


X    S    E 

0    2    0 

0    2     1 

0    2    2 

1     2    0 

1     2    2 

0     1     0 
0    2    0 
2    0    0 
2     1     0 
2    2    0 

SI 
S2 

N2 
N2 

N2 

0 
S2 

0^^ 
0^^ 
S2^ 

0 

S2 

N2 

N2 
N2 

SI 
S2 

N/A 
N/A 
N/A 

0 

S2 
N/A 
N/A 
N/A 

Explanation  of  Table  arrangement  is  provided  on  pages  8  and  9.    A  unit  on  south 
route2  blocks  only  a  unit  from  north  turning  east  (north  route2)  and  units  from  the 
south  routes  themselves.  Therefore  north  route2  has  the  priority  for  next  service  if 
there  is  no  unit  moving  on  east  route  1,  else  the  priority  will  again  be  transferred  to 
the  south  routes.    This  is  designated  by  V  in  row  (2  2  0)  and  column  (0  2  1).    Note 
that  in  Table  1,  north  route2  is  placed  ahead  of  the  south  routes  in  order  of  priority 
at  the  end  of  service  on  south  route2.    This  is  because  we  have  assumed  that  once 
service  is  completed  in  one  direction,  priority  for  next  service  will  be  given  to  the 
other  directions  which  have  been  blocked  by  the  unit  leaving  the  system.   A  unit  on 
the  east  routes,  the  last  1  or  2  of  the  triples  in  columns  2,  3,  and  5  allows  no  service 
on  south  routel  and  the  movement  on  east  routel  (1  of  column  (0  2  1))  at  the  end 
of  south  route2's  service  also  blocks  north  route2  (see  W).  The  O's  in  rows  1,  3,  and 
4  and  columns  2,  3,  and  5  are  because  units  on  east  routel, 2  continue  to  block  service 
of  units  attempting  to  move  at  the  completion  of  south  route2's  service. 
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Table  6.  SERVICE  PRIORITIES  AT  THE  END  OF  AN  EAST  UNIT  TURNING 
SOUTH  (denoted  by  east  routel):  Explanation  of  Table  arrangement  is 
provided  on  pages  8  and  9. 


X    s 

E 

0    0     1 

0    2     1 

Only  a  unit  from  the  south  moving  east 
(south  route2)  is  allowed  service  while  a 

0    0 

1 

El 

El 

0    0 

E2 

E2 

unit  from  the  east  is  moving  south  (east 

0     1 

0 

SI 

X;A 

routel).   Therefore  we  only  have  to  con- 

0    1 

1 

SI 

X/A 

sider  two  movement  combinations  just 

0     1 

2 

SI 

X  A 

prior  to  east  routel  finishing  service  and 

1     0 

0 

NT 

XT 

all  possible  combinations  of  units  wait- 

1    0 

1 

M 

XI 

ing  for  east  routel 's  service  completion. 

1     0 

2 

Nl  and  E2 

XT  and  E2 

One  or  two  units  may  receive  priority 

1    1 

0 

XI  and  SI 

X/A 

when  service  on  east  routel  is  completed 

1    1 

1 

XI  and  SI 

X  A 

according    to    the    service    policies    dis- 

1   1 

XI  and  SI 

X,A 

cussed    earlier.    The    service    on    south 

2    0 

0 

X2 

0 

route2,  the  second  2  of  (0  2  1)  in  column 

2    0 

I 

X2 

E\' 

2,  blocks  a  unit  from  the  north  waiting 

2    0 

2 

X2  and  E2 

£2^ 

to  move  east  (north  route2).    Thus  the 

2     1 

0 

X2 

x;a 

priority  initially  given  to  north  route2 

2     1 

1 

X2 

X  A 

must  be  transferred  to  the  east  routes, 

2     1 

2 

X2  and  E2 

X/A 

which  is  designated  by  V  in  this  Table. 
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Table  7.     SERVICE  PRIORITIES  AT  THE  END  OF  AN  EAST  UNIT  TURNING 
NORTH  (denoted  by  east  route2) 


N    S     E 

0    0    2 

0    2    2 

1     0    2 

1     2    2 

2    0    2 

0    0     1 
0    0    2 
0     I     0 
0     1     1 

0     1     2 

El 
E2 
SI 
SI 

SI 

El 

E2 

N/A 
N/A 

N/A 

0 

E2 
SI 
SI 
SI 

0 
E2 

N/A 
N/A 

N/A 

0^ 
E2 
0^ 
0^ 
£2^ 

Explanation  of  Table  arrangement  is  provided  on  pages  8  and  9.    A  unit  in  service 
on  east  route2  blocks  only  south  routel  and  east  routel,2.  Therefore  south  routel 
has  the  priority  for  next  service  if  there  is  no  unit  in  service  on  north  route2,  else  the 
priority  will  be  again  transferred  to  the  east  routes.  Note  that  in  Table  1,  south  routel 
places  ahead  of  the  east  routes  in  order  of  priority  at  the  end  of  service  on  south 
route2  because  we  have  assumed  that  priority  for  next  service  is  given  to  the  other 
directions  once  service  from  one  direction  is  completed.    The  unit  on  north  route2, 
the  first  2  of  column  (2  0  2)  blocks  a  unit  waiting  for  service  on  east  routel  and  south 
routel,  which  is  designated  by  V. 
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B.     ANALYTICAL  MODEL 

1.  System  Notations 

In  this  section,  we  introduce  the  notations  required  to  develop  an  analytical 
model.   Note  that  subscripts  of  «,  s,  and  e  represent  north,  south  and  east  respectively. 
north  units  Units  from  the  north  direction. 

X„  Poisson  arrival  rate  at  the  intersection  from  the  north  direction. 

fi„  Exponential  service  rate  of  a  unit  from  the  north  direction. 

Q„{t)  Number  of  north  units  waiting  for  service,  as  a  function  of  time  /. 

Next.north.route  Next  route  selected  for  the  first  unit  in  north. queue. 
north.queue  Number  of  units  waiting  for  service  from  the  north  and  next  route 

for  the  first  unit  in  the  north  queue. 
north.run  Route  that  a  north  unit  uses  now. 

P„  Probabihty  that  a  north  unit  will  take  north  route!. 

Note  that  north.queue  represents  two  quantities,  QX^)  and  next.north.route  with 
"/"  between  them.  As  an  example,  "3/1"  of  north.queue  indicates  that  there  are  3  units 
in  the  queue  and  the  first  unit  in  the  queue  is  supposed  to  use  north  route  1. 

There  are  also  south  units,  /,,  ^„  south. queue,  south.run,  next.south.route,  Q^ 
and  P,  for  the  south  direction  and  east  units,  >.„  ix„  east. queue,  east. run,  next. east. route, 
Q,  and  P,  for  the  east  direction. 

2.  A  Transition  Diagram 

Suppose  that  potential  units  arrive  in  accordance  with  Poisson  rate  /„,  ?.„  /., 
from  3  directions  and  that  service  rates  are  exponential  ^x„,  iu.„  and  n,.  We  can  then  the- 
oretically compute  the  quantities  of  interest  for  the  system;  the  mean  number  of  units 
in  the  queues  and  the  average  amount  of  time  that  an  arriving  unit  waits  in  the  queue. 
To  accompUsh  this  we  must  first  decide  upon  an  appropriate  state  space.  It  is  clear  that 
the  state  of  the  system  must  include  more  information  than  merely  the  number  of  units 
in  each  queue  of  the  system  and  the  present  movement  state  (1  of  14  movement  states 
we  discussed  earher).  For  instance,  to  decide  the  movements  after  a  service  completion, 
it  would  not  be  enough  to  specify  that  there  are  n,,  nj,  nj  units  in  each  queue  and  which 
routes  are  being  used,  but  we  would  also  have  to  know  the  routes  selected  for  the  first 
unit  in  the  queues. 

To  account  for  this,  we  define  a  state  variable  (north.queue,  south.queue, 
east.queue,  north.run,  south.run,  east. run),  where  the  first  three  entities  give  information 
about  the  number  of  units  in  the  queues  and  the  desired  routes  for  the  first  units  in  the 
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queues  and  the  last  3  entities,  north. run,  south. run  and  east. run  represent  the  routes  that 
are  presently  being  used  by  north,  south,  or  east  units.  For  example,  the  following 
sample  states  have  the  following  interpretation. 

(0/0,2/1,1/0,0,0,2): 

The  entry  0/0  means  that  there  is  no  unit  in  the  north  queue. 

The  entry  2/1  means  that  there  are  two  units  in  the  south  queue  and  south  route  1  is 
selected  for  the  first  unit  of  south  queue. 

The  entry  1/0  means  that  there  is  one  unit  in  the  east  queue.  A  route  is  not  decided 
for  the  unit  because  an  east  unit  is  still  being  served. 

The  last  three  entries  0,  0,  2  means  that  only  east  route2  out  of  the  six  routes  is  pres- 
ently being  used.  Note  it  blocks  south  route  1. 

(0/0,2/1,  1/0,  1,0,2): 

The  conditions  are  the  same  as  the  above  state  except  a  new  north  unit  uses  north 
routel  (the  0  in  position  four  has  changed  to  a  one). 

(0/0,  I/O,  1/2,  1,  1,  0)  : 

The  entry  1/0  means  that  there  is  one  unit  in  the  south  queue.  A  route  is  not  decided 
for  the  unit  because  a  south  unit  is  still  being  served. 

The  entrv'  1/2  means  that  there  is  a  unit  in  the  east  queue  and  east  route2  is  selected 
for  the  unit. 

The  entry  1  in  position  four  means  that  a  north  unit  uses  north  routel. 

The  entrv'  1  in  position  five  means  that  a  south  unit  uses  south  routel. 

It  should  be  noted  that  the  process  will  go  to  the  second  state  (0/0,  2/1,  1/0,  1, 
0,  2)  from  the  first  state  (0/0,  2/1,  1/0,  0,  0,  2)  when  a  unit  arrives  at  the  intersection  from 
the  north  and  takes  north  routel  (which  occurs  at  rate  PJ.„).  Also  the  process  will  go 
to  the  third  state  (0/0,  1/0,  1/2,  1,  1,0)  from  the  second  state  (0/0,  2/1,  1/0,  1,  0,  2)  when 
the  east  unit  completes  its  service  and  east  route2  is  selected  for  the  following  unit 
(which  occurs  at  rate  (1  —  ?,)^i,). 

As  discussed  earlier,  there  are  14  combinations  of  routes  which  are  allowed  in 
the  system.  Anyone  of  these  combinations  may  have  many  states  since  we  have  ap- 
proximately {3{n  +  1)}  states  for  each  combination  of  the  routes  where  n  is  the  length 
of  the  queues.  We  can  begin  to  analyze  the  states  from  the  following  development  of 
the  (0,0,1)  route  combination. 
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(0/0,  0/0,  0,0,  0,  0,  1) 

(0/1,  o;i,  0,1, 0,0, 1) 

(0, 2,  0/2.  0/2,  0,  0,  1) 

(1,0,1/0,  1, '0,0,0,  1) 
(1/1,  1/1,  1/1,0,0,  1) 
(1/2,  l,/2,  1/2,0,0,  1) 


(n/0,  n/0,  n/0,  0,  0,  1) 
(n/l,n/l,  n/1.0,  0,  1) 
(n'2,  n'2,  n;2,  0,  0,  1) 

Of  course  in  this  model  all  {3(«  +  1)}  states  are  not  possible  in  the  state  space 
of  the  model.  For  instance,  the  state  (1/1,  1/1,  1/1,  0,  0,  1)  is  not  in  the  state  space  since 
an  east  unit  is  still  being  served,  and  the  direction  of  the  first  unit  in  the  east  queue  has 
not  been  designated  i.e.,  the  third  entr}'  must  be  1,0  if  a  unit  from  the  east  is  being 
served.  As  a  result,  each  route  combination  may  have  less  than  {3(n  +  1)}  states  but 
this  is  still  too  many  states  to  describe  in  a  transition  diagram  for  a  whole  model.  In 
addition,  it  is  not  straightforward  to  build  a  total  state  space  since  we  would  have  to 
consider  all  interactions  of  the  routes  for  each  movement  state. 

Rather  than  a  direct  attempt  to  solve  the  whole  problem  we  can  observe  how 
the  transitions  occur  and  what  the  balance  equations  look  like  by  starting  from  the  ini- 
tial point,  that  is,  the  base  state  (0/0,  0,0,  0  0,  0,  0,  0).   A  brief  transition  diagram  that 
describes  a  small  part  of  the  whole  model  is  shown  in  Figure  2. 
3.      Solution  Procedures 

Up  to  now  we  have  decided  upon  the  state  space  and  tried  to  make  a  transition 
diagram,  although  they  are  not  fully  described.  Also  we  have  assumed  Poisson  arrivals 
and  an  exponential  service  rate  for  the  system,  so  the  model  can  be  analyzed  as  a 
continuous-time  Markov  chain,  [Ref  3:  pp.  210-273], 
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Figure  2. 


A  partial  transition  diagram  of  the  model:  The  short  arrows  at  the 
bottom  of  a  state  represent  many  more  arrows  pointing  in  and  and  out 
of  the  state.  The  notations  of  Pn,  An,  and  fin  are  used  for  p„  ,  /„,  and 
u.„.  The  notations  of  Ps,  As.  and  fis  are  used  for  the  tliree  parameters  of 
the  south  and  Pe.  /.e.  and  ij.e  are  used  for  the  three  parameters  of  the  east 
direction. 


Consider  limiting  probabiUties  pj  for  all  J  states.  Pj  =  lim  p,^{i)  where  J  is  inde- 
pendent  of  state  /  and  p^j  is  the  probability  that  the  process  enters  state  J  from  state  /  . 
If  we  assume  that  the  limiting  probabilities  Pj  exist,  then  they  can  be  obtained  by  the 
following  formula,  [Ref.  4:  pp.  246-248]. 
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^j Pj  =  ZJk  Pkj Pk    for  all  state  j,  k.  (2.1) 

/Pj=\,     for  ail  State  J.  where  (2.2) 

J 

Vjpj  :  Rate  at  which  the  process  leaves  state  y. 

y^kPkjPk  '•  ^ate  at  which  the  process  enters  state  J. 

V,-  :  Leaving  rate  from  state  y. 

Vj^  :  Leaving  rate  from  state  k. 

Pj  :  Proportion  of  time  the  process  is  in  state  y. 

Pi^  :  Proportion  of  time  the  process  is  in  state  k. 

Pi^j  :  Probabihty  of  transitions  from  state  k  to  state/ 

Once  we  solve  the  balance  equations,  it  is  easy  to  answer  our  questions  in  terms 
of  these  limiting  probabilities.  However,  as  we  see  in  Figure  2,  the  number  of  states  de- 
rived from  the  base  state  increases  geometrically.  Suppose  the  maximum  number  of  "n" 
is  9,  then  we  have  approximately  (10  x  3)  x  14  =  415,800  states  and  the  same  number 
of  balance  equations  where  14  indicates  the  number  of  all  possible  combinations  of 
movements.  It  is  not  straightforward  to  develop  the  transition  diagram  and  to  write  all 
the  balance  equations  and  as  a  result,  the  approach  to  solving  the  problem  necessitates 
the  simulation  method. 
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III.      SIMULATION  MODEL  1 

A.  ASSUMPTIONS  AND  NOTATIONS 

The  assumptions  and  notations  of  simulation  iModel  1  remain  the  same  as  those  for 
the  analytical  model  discussed  in  Chapter  II,  except  that  the  interarrival  and  service 
times  may  now  have  a  Gamma  distribution  with  various  parameters,  [Ref  5:  pp. 
267-273]. 

B.  INPUTS  TO  THE  MODEL 

Our  primary  goal  is  to  examine  the  waiting  time  and  the  number  of  units  in  the 
queue  of  the  steady  state  under  certain  conditions  given  by  specified  input  parameters. 
The  simulation  model  is  made  ready  to  run  by  entering  the  necessary  input  values.  This 
part  of  the  model  is  extremely  important  since  it  provides  the  user  an  opportunity  to 
specify  the  particular  conditions  which  are  to  be  examined,  as  well  as  the  values  which 
establish  the  boundary  conditions  of  the  simulation.  These  input  values  represent  fea- 
tures of  the  system  such  as  the  probability  of  taking  route!  and  the  shape  of  the  service 
and  interarrival  time  distributions.  Hence,  for  Model  1,  we  have  a  total  of  three  factors 
which  we  divide  into  different  levels  in  Table  8  for  the  purpose  of  controlling  the  size  of 
the  simulation. 

We  assume  that  successive  service  and  interarrival  time  distributions  for  the  three 
directions  are  independent  and  identically  distributed.  In  addition,  we  assume  the 
probability  of  taking  route  1  is  the  same  for  each  direction.  Thus  from  Table  8  we  have 
2  levels  of6^,k^  and  ki,  3  levels  of  P,  and  3  levels  of  traffic  intensity.  Thus  we  initially 
simulated  (levels  of  ^J  x  (levels  of /cj  x  (levels  of /cj)  x  (levels  of  P,)  x  (levels  of  p)  = 
2x2x2x3x3  =  72  cases  of  svstem  conditions.  Note  that  we  have  4  different  com- 
binations  of  the  Gamma  shape  parameters  /cj  and  /cj;  (1,1),  (1,4),  (4,1),  (4,4).  The  value 
1  for  /:,  or  /cj  represents  the  case  where  service  time  or  interarrival  times  are  skewed,  i.e. 
the  arrivals  form  a  Poisson  process  and  the  service  rate  is  exponential,  as  in  the  analyt- 
ical model.  The  value  4  for  the  shape  parameter  indicates  that  the  Gamma  distribution 
is  regular,  [Ref  6:  p.  VIII. 7]. 
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Table  8.     MAJOR  FACTORS  AND  LEVELS  OF  SIMULATION. 


Factors 

Levels 

Service  time  distributions 

(denoted  by  the  subscript  1 

on  parameters  6  and  k) 

Gamma(0i,/c,)  where 
0,  =  2.5,  5.0  and  A,  =  1,4 

Trafiic  intensities  p 
r     _       mean  service  time     ^ 

p  =  0.5,  0.4,  0.33 

^^        mean  interarrival  time  ' 

Interarrival  time  distributions 

(denoted  by  the  subscript  2 

on  parameters  6  and  A) 

Gamma(^2,A2)  where 
6,  =  ^  and  ^2=  1,4 

Probabilities  of  taking  route  1 

Pi  =  0.25,  0.5,  0.75 

C.      MODEL  IMPLEMENTATION 

\.      Simulation  Language 

The  model  is  written  in  Simscript  11.5,  which  has  become  the  standard  simu- 
lation language  for  miUtar}'  apphcations  and  also  a  powerful  and  flexible  simulation  tool 
in  all  areas  of  miUtar>'  modeling.  Simscript  II. 5  reduces  simulation  programming  time 
and  cost  several  fold  compared  to  FORTRAN,  and  is  enormously  more  efficient  and 
flexible  than  interpretive  simulation  techniques,  [Ref  7,81. 
2.  Model  Structure 
a.      General 

The  program  contains  a  preamble,  a  main,  4  routines  and  9  event  routines. 
The  4  routines  include  read.data,  initialize,  single.run. report,  and  sum.report.  The  9 
event  routines  are  north(south, east). unit. generator,  ne\v.north(south,east).unit,  and 
north(south,east). service.  Routines  are  called  by  the  main  routine  or  by  other  routines 
and  return  control  to  the  calling  routine.  Event  routines  are  special  types  of  subroutines 
which  are  accessed  by  being  scheduled,  [Ref  9,10].  The  functions  of  each  routine  and 
interactions  among  routines  are  as  follows: 
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b.  Preamble 

It  contains  the  definitions  of  entities,  attributes  and  sets,  and  other  global 
information. 

c.  Main 

This  is  the  first  routine  to  receive  control  during  the  execution  of  a  simu- 
lation run.  Its  function  is  to  initialize  the  system  and  to  combine  the  input  parameters 
for  performing  a  series  of  simulations.  It  also  computes  and  prints  some  simulation  sta- 
tistics. 

d.  Read.data 

This  routine  reads  all  the  input  parameters  provided  by  a  user.  A  user 
specifies  these  parameters  in  the  main  routine.  The  input  parameters  are  three  seed 
stream  numbers,  the  desired  number  of  units  to  be  simulated,  the  desired  number  of 
replications,  the  various  Gamma  distribution  parameters  for  the  interarrival  and  service 
times,  and  the  probability  of  taking  route!  for  each  direction  of  the  intersection.  In  this 
thesis  we  examine  only  the  cases  described  in  Table  8. 

e.  Initialize 

The  spaces  for  all  array  elements  are  allocated.  It  also  schedules  initial 
events. 

/.      Single. run.report 

All  information  of  a  single  simulation  run  is  collected  here  and  a  verification 
of  the  model  can  be  doi'e  by  examining  this  report.  A  sample  report  is  shown  in  Table 
11  of  Appendix  A. 

g.      Sum.report 

This  routine  summarizes  the  output  from  single. run. report.  The  report  can 
be  used  to  note  the  stationary'  statistics  and  to  observe  how  the  output  statistics  change 
with  respect  to  the  input  parameters.  In  addition  it  can  be  used  to  analyze  the  boundary 
conditions  for  stability;  this  being  one  of  the  research  objectives.  The  study  of  this  re- 
port is  used  to  determine  an  appropriate  design  size  for  Model  1.  A  sample  section  of 
the  resulting  report  is  shown  in  Table  12  of  Appendix  A. 
h.      North(  south,east )  .unit. generator 

Note  that  south  and  east  in  the  parenthesis  indicate  that  the  same  routines 
also  exist  for  the  south  and  east  directions.  New.north.unit  event  routine  is  scheduled 
immediately  after  the  north. unit. generator  event  routine  is  called  and  this  event  itself  is 
rescheduled  for  a  next  arrival.     The  interarrival  time  for  the  next  north  arrival  is  a 
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Gamma  random  variate  whose  parameters,  as  described  previously,  are  inputs  to  the 
model.  The  termination  of  simulation  is  also  controlled  by  this  event  routine  after  the 
desired  number  of  units  have  been  simulated  from  each  of  the  three  directions. 
I.      New.north(  south,east )  .unit 

A  north  unit  is  created  immediately  after  this  event  routine  is  called  and  this 
unit  is  either  served  directly  or  filed  in  the  queue.  Whether  it  is  served  or  filed  in  the 
queue  and  which  route  is  selected  for  the  unit,  depends  on  other  movements  at  the 
intersection  and  the  route  selection  procedures  we  discussed  in  Chapter  II  (see  route 
selection  procedure  on  page  4).  If  the  unit  is  served  immediately,  a  north. service  event 
routine  is  scheduled  to  occur  at  the  end  of  service.  Note  that  the  service  time  is  a 
Gamma  random  variate. 

j.      North(  south, east). service 

A  unit  that  has  passed  the  intersection  is  destroyed  immediately  and  this 
unit  leaves  the  system.  The  units  blocked  by  the  leaving  unit  are  served  according  to  the 
service  priorities  discussed  in  Chapter  II  and  the  appropriate  north(south, east). service 
event  routines  will  be  scheduled  based  on  new  service  times  (Gamma  random  variates). 
3.     Logic  of  the  Model 

The  flow  charts  of  4  major  event  routines  for  the  north  and  a  macro  flow  chart 
for  the  whole  model  are  depicted  in  Figure  3  through  Figure  7  respectively.  The  flow 
charts  for  other  directions  are  omitted  to  avoid  reduiidancy.  The  logic  of  the  event  rou- 
tines for  the  south  and  east  is  exactly  the  same  as  that  for  the  north,  though  the  condi- 
tions for  movement  to  a  given  route  is  quite  different. 

We  can  divide  the  north(south,east). service  event  routine  into  two  phases;  a 
route  selection  phase  for  the  first  unit  in  the  queue  and  a  decision  phase  for  the  next 
movements.  There  is  no  diflerence  between  the  north,  south  or  east  directions  in  the 
route  selection  phase.  However,  the  decision  phase  is  quite  diflerent  for  each  direction 
since  each  route  blocks  diflerent  routes  and  the  conditions  of  movement  to  a  blocked 
route  is  diflerent.  As  should  be  apparent  at  this  point,  the  decision  phase  of  the 
north(south,east). service  event  routine  must  correspond  with  the  service  pohcies  and 
priority  rules  discussed  earlier.   The  program  hsting  is  presented  in  Appendix  B. 
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Figure  3. 


Flow  chart  for  north.unit.generator  event  routine:  Jn  is  the  number  of 
north. units  that  have  passed  the  intersection,  js  is  for  south. units  and  je 
is  for  east. units.  NI  is  the  minimum  desired  number  of  units  to  be  sim- 
ulated from  each  direction  and  is  selected  by  a  user.  This  routine  will  be 
rescheduled  as  long  as  at  least  one  of  jn,  js  or  je  is  less  than  Nl  and  will 
stop  being  rescheduled  if  all  three  values  are  greater  than  or  equal  to 
Nl. 
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Figure  4. 


Flow  chart  for  new.nortli.unit  event  routine:  U  is  a  Uniform  (0,1)  ran- 
dom variate  and  is  compared  wirh  P,  for  route  selection  of  tlic  unit  where 
P,  is  the  probability  of  taking  route  I.  An  arriving  unit  which  is  created 
in  this  routine  can  be  served  immediately  or  filed  in  the  queue  depending 
upon  queue  condition,  route  selection  and  the  other  movements  at  the 
intersection. 
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Figure  5.  Flow  chart  for  nortii.service  event  routine:  U  is  a  Uniform  (0,1)  random 
variate  and  Pj  is  the  probability  of  taking  route  1.  In  this  routine,  prior- 
ities for  service  after  a  service  completion  on  north  route  1  or  2  are  given 
to  appropriate  routes  according  to  the  service  policies  discussed  in  Table 
2  through  Table  7  in  Chapter  II. 


Figure  6. 


Flow  chart  for  nortli.service  event  routine  (continued):  The  priorities  for 
next  service  after  service  completion  in  the  south  or  east  are  decided  in 
the  south(east). service  routines.  Note  that  in  the  program  the  value  of 
next. north(  south. east). route  is  changed  to  0  whenever  a 
north(south. east). service  event  routine  is  scheduled. 
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Figure  7. 


A  macro  flow  chart  for  the  model:  \  is  the  smallest  number  of 
north(south,east)  units  that  have  passed  the  intersection  and  Nl  (se- 
lected by  the  user)  is  the  desired  number  of  units  to  be  simulated.  R  and 
NR  indicate  the  number  of  replications  completed  and  the  desired  num- 
ber of  replications  respectively.  C  and  NC  are  used  to  check  if  all  simu- 
lation conditions  are  examined  where  C  indicates  the  number  of 
conditions  that  have  been  examined,  and  NC  is  the  total  number  of 
conditions  to  be  examined.  See  Table  S  for  the  simulation  conditions. 
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D.      RESULTS 

1.  Statistical  Considerations 

When  the  combinations  of  many  different  input  specifications  to  a  simulation 
are  to  be  evaluated,  as  in  this  problem,  we  must  first  decide  upon  which  combinations 
of  factors  we  wish  to  spend  our  limited  simulation  resources.  So,  many  preliminar\'  ex- 
periments using  various  combinations  of  input  parameters  were  investigated,  and  those 
conditions  under  which  the  system  was  not  stable  were  used  as  boundary  conditions  in 
the  experimental  design.  As  a  result,  a  finite  number  of  discrete  level  experiments  were 
conducted  as  a  series  of  essentially  independent  tests  in  looking  for  trends  in  each  di- 
rection. We  used  100  replications  within  each  run  to  achieve  sufficient  precision  for  each 
response  estimate. 

In  order  to  reduce  the  variance  among  simulation  runs,  common  random  num- 
ber sequences  were  used  for  all  the  simulations  of  different  factor  combinations,  [Ref. 
11:  pp. 250-252].  Three  separate  sequences  were  used  for  each  simulation  (interarrival 
time,  service  time,  and  route  selection).  The  random  number  sequences  were  continued 
to  the  next  replication,  so  that  there  was  no  overlap  of  the  random  number  sequences 
and  we  could  therefore  assume  the  replications  for  each  set  of  factors  were  independent 
of  each  other. 

2.  Simulation  Outputs 

Whenever  we  have  one  replication  of  a  simulation,  as  discussed  before,  a  par- 
ticular realization  of  the  output  sequences  will  be  generated.  Examples  of  random  se- 
quences of  interest  would  be  a  running  average  of  successive  waiting  times  and  a  running 
average  of  queue  length  experienced  by  the  first  r  north( south, east)  units  for  r  = 
1,...,N1,  where  NT  is  the  desired  number  of  units  to  be  simulated  from  each  direction. 
Nl  is  a  user  input  and  must  be  large  enough  to  provide  stability. 

For  this  study,  the  following  sequences  of  interest  were  collected  and  computed 
in  the  single. run.report  and  sum.report  routines. 

Wi^,    :    Waiting  time  for  the  r""  departing  unit  from  the  /"'  direction,  where  /  can  have 
values  north,  south  or  east. 

Ti,.      :    Departure  time  for  the  r"'  unit  from  the  /"'  direction. 

Qi,      :    The  number  of  units  waiting  in  the  i"'  direction  at  the  /-"'  unit's  departure. 

A  Wi^,   :    The  running  average  of  the  first  r  waiting  times  obtained  by  the  formula: 
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r 

AWi,  =  -^2_j^^i,v^      /=north,south,east  andr=  1,...,N1.  (3.1) 

v=l 

AQi^   :    The  running  average  of  the  number  of  units  waiting  in  the  /"'  direction  during 

the  time  0  <  /  <  T,,.  It  is  a  function  of  time  and  a  random  sequence,  but  it  can  also  be 
viewed  as  a  sample  of  random  processes  defined  on  continuous  time  or  as  a  function  or 
integral  of  such  a  process.  For  example,  if  we  let  Q,{t)  be  the  number  of  units  waiting  in 
the  i"'  direction's  queue  at  time  i,  then  {Q,{t);  /  >  0}  is  a  random  process  with  a  contin- 
uous parameter  /.  For  /  >  0,  AQ,^  can  be  expressed  by  the  following  formula,  [Ref.  11: 
p.272-273]: 


T 


0-^IJ.)diJ.,       /  =  north,south,east  and  r  =  1,...,N1.  (3.2) 


which  is  easily  cornputed  by  the  ACCUMULATE  function  in  Simscript  II. 5,  [Ref.  7: 
p.272-273]. 

The  sequences  of  W^,  ,  A  Jr,_,  ,AQ,^,  and  the  process  Q,{t)  for  the  north  direction 
and  100  units  of  Nl  with  skewed  interarrival  and  service  time  distributions,  traffic 
intensity  =  0.5  and  Pi  =  0.75  are  shown  in  Figure  8.  Note  that  this  set  of  simulation 
conditions  is  also  used  for  further  system  output  displays. 
3.      Characteristics  of  Output  Sequences 

As  we  discussed  above,  the  outputs  of  a  simulation  are  sequences  of  random 
variables.  In  this  study  we  are  interested  in  certain  probability  distribution  characteristics 
and  our  objective  is  to  estimate  them  via  this  simulation  since  they  are  unknown  and  can 
not  be  computed  with  theoretical  analysis. 

Assume  we  have  M  independent  observations  of  the  sequence  [A  W,J  and 
{AQ^,},  i=  north, south,east  and  r  =  1,...,N1  through  independent  replications  of  a  sim- 
ulation. To  obtain  independent  replications  of  a  simulation,  we  generate  independent 
sequences  of  random  variables  to  drive  the  simulation.  These  sequences  are  recursively 
generated,  as  discussed  previously  in  this  section. 

To  estimate  the  mean  values  and  obtain  confidence  intervals  of  the  observations 
from  sets  of  sample  sequences  (such  as  depicted  in  Figures  9  and  10)  we  let  ^x^^^.  and 
ol^  be  the  mean  and  variance  sequences  of /4H^,,  and /i^o  and  ct^^.  be  those  of /i 2, „ 
i.e.. 
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Figure  8.      Output  sequences  obtained  from  a  single  simulation  run 


^AW^.^^^'^^^lrl 


^A0„=  ^l^Qhr'l^ 


i  =  north, south, east  and  r  =  1,...,N1. 


/  =  north,south,east  and  r  =  1,...,N1. 


1  '' 

a        =  EKAlVii-  —  n       )"],        /  =  north,south,east  and  r=  1,.,.,N1, 


1  1 

a  .^   =  E[{AQ,.  — fx ^^)'\,        /=  north,south,east  and  r  =  1,...,N1, 


(3.3) 
(3.4) 
(3.5) 
(3.6) 


^^^  =  A  Wi,.  =  —r-r  y^A iVi^r,m  »      ^  =  north,south,east  and  r  =  1,...,N1.       (3.7) 
Thus,  the  sequence  of  standard  point  estimators  of  fx        is  the  sequence  of 


sample  means  ju^^    ;  (3.7)  where  M  is  the  desired  number  of  replications  of  the  simu- 
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lations  and  AW^,„  is  the  value  of  the  m'*  replication  of  AW,,.  A  confidence  interval  for 
the  sequence  ij..^    is  provided  with  the  t-statistic  as  follows,  [Ref  11:  pp.282-283]. 

M 

Sa  w„  =  S,      =  -^  YJ-^  W„_„  -  AW  J,  (3.8) 

i  =  north,south,east  and  r=  1,...N1  and  M  =  number  of  replications 


A 


5.„.   =o^...       /=  north.south.east  and  r  =  1,...,N1.  (3.9) 

A  iV,,.  A  fy,, 

Var{jji       ]  =  — 7T~  ,     /=  north,south,east  and  r=  1,...,N1.  (3.10) 

Given  the  above  definitions  a  random  variable, 

A 

^/l  W, ,  ~  ^A  IV, , 


^A  »;, 
has  approximately  a  t-distribution  with  M  —  1  degrees  of  freedom  and  we  have 

A 

Prob{t^^_,ial2)  <  ^^"^-^     ^^"^--  <  r,,_,(l  -  a/2)}  =  1  -  a,  (3.11) 

'^A  »;, 

where  /=  north, south, east,  r=  1,...N1,  M=  number  of  replications,  and  /M-i(a/2)  is  the 
a/2  quantile  of  the  t-distribution  with  M  —  1  degrees  of  freedom.  Since  the  t-distribution 
is  symmetric  about  zero,  we  obtain, 

S  S 

Pr«*{;i^^^  -,„.,(!  -a/2)-^<^^^  <M,^,^  +  r„_,(I  -a/2)-^}  (3.12) 

«  1  —  a,    /  =  north,south,east  and  r  =  1,...,N1. 


S  5 

A  AW  A  W 

Therefore  the  two  values,  ix^^^^  -  t^_S\  -  a/2)  —^  and  li^^^^  +  rw_,(l  -  a/2)  -^  are 
100(1  -  a)%  confidence  intervals  for  /x       .  Given  the  previously  discussed  parameter 
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conditions,  the  realizations  o[  li^^  and  the  confidence  intervals  of  u^,^  obtained  from 
the  6  sequences  in  Figure  9  and  Figure  10  are  depicted  in  Figure  11  and  12  for  the  north 
and  south  directions  respectively. 

Let  u,„,      be  the  500"'  value  of  the  sequence  u,^  ;  the  mean  estimator  of  the 
running  average  of  successive  waiting  times  experienced  by  the  first  500  units  and  let 


AW, 


i\500 


be  the  standard  deviation  estimator  of  a 


AV, 


1.500 


The  values  of  pt 


AW, 


1.500 


AW 


1.500 


ob- 


tained from  various  simulation  conditions  are  shown  in  Tables  16  through  18  of  Ap- 
pendix A.  The  values  o'l  Ix^^,  for  different  parameter  sets  and  different  directions  will 
be  analyzed  in  detail  in  Chapter  V, 


Figure  9. 


Six  realizations  of  an  output  sequence  (north  queue):  This  figure  shows 
six  realizations  of  the  sequence  {A  /F,,}  for  value  i  corresponding  to  the 
north  direction.  We  can  see  the  random  appearance  and  the  variation 
from  sequence  to  sequence.  The  distribution  characteristics  of  the  ran- 
dom sequences  such  as  these  are  what  we  are  interested  in  estimating. 
These  characteristics  are  estimated  from  the  sequences  in  the  figure  and 
are  shown  in  Figure  11. 
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Figure  10.  SLx  realizations  of  an  output  sequence  (south  queue):  This  figure  shows 
six  realizations  of  the  sequence  [AWi^,]  for  value  i  corresponding  to  the 
south  direction.  The  sequences  in  this  figure  also  show  random  ap- 
pearance and  variation  from  sequence  to  sequence.  The  distribution 
characteristics  of  the  random  sequences  are  estimated  and  shown  in 
Figure  12. 

4.      Stability  and  Convergence  Considerations. 

Note  that  there  are  three  separate  queues  in  the  system.  As  we  see  in  Figure 
11  and  Figure  12,  the  sequence  [li^^  }  seems  to  converge  to  n.^  in  the  north  and  south 
queues  under  the  given  traffic  conditions.  We  found  out  that  the  east  queue  is  not  stable 
under  the  same  conditions.  We  know  that  it  makes  a  difference  in  simulation  difficulty 
if  the  sequence  of  interest  settles  down  and  converges  to  a  stationar>'  distribution  after 
an  initial  transient  phase.  If  the  sequence  settles  down,  then  we  may  only  need  to  con- 
sider the  distribution  that  exists  after  the  transient  phase  and  not  a  separate  distribution 
for  each  value  in  the  sequence,  [Ref  6:  p. II  1.22]. 

To  test  convergence  to  the  stationary'  distribution,  we  made  preliminarv',  inde- 
pendent replications  of  the  simulation  over  a  range  of  r  (number  of  units  simulated) 
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Mean  and  C.I  sequences  of  waiting  time  (north  queue):  This  figure 
shows  a  realization  of  [ix.,..  )  and  confidence  intervals  of  (ju.„,  }  for 
value  i  corresponding  to  the  North  direction.  The  sequence  {ii^^,  ]  is 
a  random  sequence  obtained  from  6  sequences  o^  AW',,  in  Figure  9. 
Since  E\Ji^^,  ]  =M^^,  ,  it  fluctuates  about  the  converging  sequence  of 
means.  The  larger  the  value  of  M  the  narrower  the  fluctuations  and  the 
closer  the  approximation  to  the  sequence  [ii.^  ].  Hence,  by  taking  a 
large  value  M  (100  in  this  thesis)  a  reasonable  estimate  can  be  made 
on  a  convergence  point  for  the  sequence  {/j.^^  }.  Notice  that  from  this 
figure  a  reasonable  judgement  can  be  made  as  to  the  extent  of  the 
transient  phase,  whereas  this  would  be  impossible  from  any  individual 
output  sequence.  In  this  figure  a  reasonable  estimate  for  the  end  of  the 
transient  phase  appears  to  be  about  the  first  100  waiting  times. 


sufiicient  to  observe  convergence.  The  sequence  {/j_,  }  was  estimated  by  generating  a 
sequence  of  sample  means  as  discussed  before.  The  most  straightforward  convergence 
test  procedure  would  be  to  collect  and  plot  the  whole  sequence  and  check  the  conver- 
gence by  examining  the  trends.  In  practice,  however,  such  an  investigation  requires  more 
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Figure  12.  Mean  and  C.I  sequences  of  waiting  time  (south  queue):  This  figure 
shows  a  realization  of  {li^^  }  and  confidence  intervals  of  {/j..^  }  for 
value  i  corresponding  to  the  south  direction.  The  sequence  {/i^^,  }  is  a 
random  sequence  obtained  from  6  sequences  of  AW^,  in  Figure  10.  In 
this  figure  we  can  see  that  the  south  queue  is  stable  and  a  reasonable 
estimate  for  the  end  of  the  transient  phase  appears  to  be  about  the  first 
300  waiting  times. 

computing  time  and  effort  than  would  be  reasonable.  Therefore  w^e  collected  ever\'  lO' 
data  point  after  the  first  300  points  of  a  sample  mean  sequence,  such  as  the  one  seen  in 
Figure  11,  to  save  our  limited  simulation  resources.  To  examine  convergence  of  the  re- 
sulting sequence  we  manually  compared  the  values  over  the  sequence  range. 

The  convergence  tests  revealed  the  system  is  partially  stable  in  the  north  and 
south  queues  for  Model  1  under  the  condition  that  traffic  intensity  =  0.5,  P,  =  0.75  and 
the  ratio  of  mean  interarrival  time  for  east  over  mean  interarrival  time  for  north  or  south 
is  1.  To  illustrate  the  east  direction's  unstable  queue,  we  collected  6  complete  sequences 
of /4W^,„,,,  and  plotted  them  in  Figure  13.  Note  the  upward  trend  in  the  sequences.  As 
we  increase  the  ratio;  (mean  interarrival  time  for  east)  /  (mean  interarrival  time  for 
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north),  from  I  to  1.3  by  increments  of  0.1,  the  east  queue  becomes  more  stable  (Figure 
14). 


Figure  13.  Six  realizations  of  an  output  sequence  (east  que'ie):  This  figure  shows 
six  reahzations  of  the  sequence  [A  U\^]  for  value  i  corresponding  to  the 
east  direction.  The  sequences  show  an  upward  trend  and  indicate  there 
is  no  steady  state  distribution  for  the  mean  of  the  sequences.  The 
characteristics  of  each  data  point  in  tlie  mean  sequence  can  be  esti- 
mated from  the  corresponding  data  points  in  the  six  A  W^^,  sequences. 

5.      Data  Collection  and  Analysis  Tools 

We  have  powerful  tools  for  data  collection  and  analysis  such  as  "print"  and 
"accumulate/tally"  statements  in  PC-simscript.  These  statements  can  be  used  for  data 
collection  and  some  statistical  calculations.  Other  software  packages  such  as 
STATGRAF  in  PC  can  be  used  for  more  sophisticated  data  analysis.  In  addition  to 
these  tools  for  use  on  the  personal  computer,  we  used  IBM  303 3 AP  mainframe  support 
to  efficiently  manipulate  much  of  the  data.  Data  was  collected  on  floppy  diskettes  and 
transferred  to  the  IBM  mainframe  with  the  3270  PC's  Host-Supported  File  Transfer 
program.    Once  transferred  to  the  mainframe,  we  edited  the  data  file  prior  to  analyzing 
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the  data  with  GRAFSTAT.  GRAFSTAT  graphical  displays  were  extremely  useful  for 
checking  the  characteristics  of  large  data  sets  such  as  required  for  Figure  9  and  Figures 
16  through  22  in  Chapter  V. 
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Figure  14. 


Stability  test  on  the  east  queue:  We  tested  the  stability  of  the  east 
queue  by  increasing  the  ratio  of  mean  interarrival  time  for  the  east  over 
the  mean  interarrival  time  for  the  north  from  1.0  to  1.3.  The  estimators 
of  E\.AJV,_,']  =  {n^^  )  is  ll^^^,  ,  the  sequence  of  the  sample  means  of 
{A  fF,  J  for  value  i  corresponding  to  the  east  direction.  The  east  queue 
becomes  more  stable  as  the  "ratio"  increases. 
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IV.      SIMULATION  MODEL  2 

A.      ASSUMPTIONS  AND  CHANGES  TO  THE  MODEL 

After  developing  a  model  of  the  system  using  the  described  assumptions,  we  propose 
and  simulate  a  road  structure  to  reduce  the  restrictions  of  the  current  road  system. 
Figure  15  depicts  the  proposed  system  which  will  be  addressed  as  Model  2.  It  is  obvious 
that  Model  2  allows  more  combinations  of  routes,  so  that  traffic  performance  at  the 
intersection  should  be  increased.  Actually  4  combinations  of  routes,  (0,1,1),  (2,0,1), 
(2,1,0),  and  (2,1,1),  are  added  to  the  14  cases  of  Model  1.  In  case  of  blockages,  the  pri- 
orities for  next  service  arc  shown  in  Table  9.  Note  that  in  Model  2,  each  route  only 
blocks  the  route  that  merges  with  it  in  a  common  direction.  All  other  assumptions  of 
Model  2  remain  the  same  as  in  Model  1. 

Table  9.  PRIORITIES  AT  THE  END  OF  SERVICE  OF  A  BLOCKING  ROUTE 
(MODEL  2):  See  Table  1  and  explanation  of  Table  arrangement  pro- 
vided on  page  7  for  comparison  with  Model  1. 


Route  which  is  about 
to  finish  service 

Next  Priority,  in  order,  at  the  end  of  blocking 

north  route  1 
north  route2 
south  route  1 
south  route2 
east  route  1 
east  route2 

east  routel,  (north  routel,  north  route2) 
south  route2,  (north  routel,  north  route2) 
east  route2,  (south  routel,  south  route2) 
north  route2,  (south  routel,  south  route2) 
north  routel,  (east  routel,  east  route2) 
south  routel,  (east  routel,  east  route2) 
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Figure  15. 


A  proposed  modification  of  the  current  system  (Model  2):      Here  several 

blockages  have  been  removed.  In  particular  tralfic  coming  from  east 
routel  is  never  blocked  by  traffic  coming  from  south  routel  though  it 
still  has  to  merge  with  traffic  from  north  routel.  Similarly,  traffic  on 
north  route2  is  no  longer  blocked  bv  either  traffic  on  east  routel  or 
south  routel,  though  it  still  can  be  blocked  by  traffiic  from  South 
route2.  Finally,  traffiic  on  south  routel  is  no  longer  blocked  by  traffiic 
on  north  route2  or  east  routel  but  can  still  be  blocked  by  traffiic  from 
east  route2. 
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B.  INPUTS  TO  THE  MODEL 

Inputs  are  the  same  as  in  Model  1. 

C.  MODEL  IMPLEMENTATION 

Since  the  only  difierence  between  Model  1  and  Model  2  are  the  additional  route 
combinations,  all  other  events  and  routines  remain  the  same  as  the  first  model  except  the 
events  related  to  evoking  service.  These  routines  are  new.north(south,east).unit  and 
north(south,east). service  event  routines.  As  an  example,  consider  the  route  combination 
(2,1,0);  it  is  not  allowed  in  Model  1  because  of  the  movements  of  north  route2  and  south 
route  1  are  in  conflict,  however,  the  two  movements  are  allowed  in  Model  2.  Thus  the 
decision  phase  of  the  event  routines  are  revised  for  the  new  model  as  required  to  allow 
the  additional  route  combinations.  The  program  listing  is  shown  in  Appendix  C. 

D.  MODEL  VERIFICATION  AND  VALIDATION 

Having  set  up  the  service  policies,  system  assumptions  and  developed  the  simulation 
model,  the  next  step  is  to  verify  the  model;  that  is,  to  detect  and  correct  discrepancies 
between  the  intended  and  actual  execution  processing  performed  by  the  model. 

A  simulation  model  of  a  complicated  process  such  as  this  model  can  not  be  de- 
bugged or  verified  unless  comprehensive  diagnostic  outputs,  which  completely  describe 
the  operation  of  the  simulated  system,  are  provided.  Diagnostics  are  needed  only  for 
debugging,  being  much  too  cumbersome  and  costly  to  use  in  full-scale  runs.  The  most 
powerful  and  economic  diagnostic  tool  used  in  this  model  is  output  which  describes  each 
simulated  event  on  a  separate  line  with  almost  all  relevant  information  about  each  unit's 
status.  The  data  collected  in  the  single. run.report  routine  is  also  used  for  verification 
of  the  model.  The  sample  sections  of  this  diagnostic  and  the  data  collected  in  the 
single. run. report  routine  are  shown  in  Tables  13  through  15  of  Appendix  A.  Another 
verification  tool  used  was  to  examine  the  logic  of  the  flow  charts  discussed  earlier. 

To  validate  the  model  as  a  sufficiently  accurate  representation  of  the  intersection 
operations  for  the  purpose  of  the  study,  that  is,  to  answer  the  question  of  whether  there 
is  a  correspondence  between  the  model  and  the  physical  intersection  system,  would  re- 
quire the  models  results  be  compared  with  that  of  the  real  system  in  a  given  situation. 
It  would  be  impossible  to  create  the  physical  system  conditions  simulated  by  the  model 
because  of  the  many  restricted  military  situations  in  R.O.K.  We  can  not  exercise  the 
movements  of  all  forces  to  simulate  wartime  traffic  given  the  current  situation.  In  spite 
of  the  restrictions,  manual  checks  of  the  model  flow  charts  and  careful  examination  of 
the  assumptions  and  service  policies  may  be  conducted  to  confirm  that  the  model  cor- 
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responds  to  the  actual  system.  If  we  assume  that  military  police  comrol  the  traffic  at  the 
three  way  intersection,  the  general  assumptions  and  service  poUcies  in  the  models  do  not 
correspond  to  the  real  world,  however,  that  is  not  a  problem  since  we  are  interested  in 
overall  traffic  performance  at  the  intersection  and  service  priority  issues  do  not  affect 
our  study  objective. 
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V.      ANALYSIS 

A.      SIMULATION  RESULTS 

1.      Current  Road  System  (Model  1) 

We  present  a  graphical  approach  in  organizing  and  presenting  the  simulation 
results,  [Ref.  6:  p.  VIII. 3-23].  The  interest  is  showing,  in  a  simple  graphical  form,  how 
the  chosen  quantification,  i.e.  mean  waiting  times  at  the  intersection,  var\-  with  factors 
such  as  traffic  intensities,  service(unit)  length  and  interarrival  time  distributions.  Figures 
16  through  18  show  the  results  of  the  first  model  with  factor  Pi  =  0.75,  0.50  and  0.25 
respectively. 

In  Figure  16,  we  have  six  sub-figures;  three  in  the  upper  half  with  mean  service 
times  of  2.5  for  all  directions  and  three  in  the  lower  half  with  mean  service  times  of  5.0 
for  all  directions.  Each  cell  in  the  sub-figure  is  indexed  by  the  traffic  intensity  and  the 
combinations  of  interarrival  and  service  time  distributions  defined  here  to  be  SS,  SR, 
RS,  and  RR,  (see  Figure  16  for  an  explanation  of  these  terms)  and  contains  the  esti- 
mated value  of  Ai,„,  (see  Equations  3.7  and  3.9).  In  Figures  17  and  18,  we  can  observe 
mean  waiting  times  of  the  three  directions  for  all  combinations  of  interarrival  and  service 
time  distributions,  given  traffic  intensity  =  0.5  and  mean  service  time=  2.5.  In  addition. 
Figure  22  shows  the  changes  to  the  mean  waiting  times  of  the  east  queue  together  with 
changes  to  the  mean  interarrival  time  ratios  for  the  three  directions,  given  Model  1  with 
P,  =  0.75,  traffic  intensity  =  0.5,  and  mean  service  time=  2.5.  Note  that  the  waiting  times 
in  Figures  17  and  18  should  be  compared  with  the  waiting  times  on  the  graphs  in  Figure 
16  with  mean  service  time  2.5  and  traffic  intensity  of  0.5.  Note  too  that  the  vertical 
scales  in  Figures  16  through  18  are  different.  It  would  be  preferable  to  use  the  same  scale 
to  give  a  depiction  of  the  relative  difference  of  the  mean  waiting  time  as  simulation  pa- 
rameters change,  however,  this  would  require  either  a  very  large  display  or  a  transfor- 
mation of  the  data  (such  as  a  log  scaling).  Neither  of  these  were  desirable,  thus  we  need 
to  look  at  the  vertical  scales  when  comparing  the  estimators  between  sub-figures. 

Hence,  it  is  relatively  easy  to  visually  inspect  Figure  16  to  determine  the  effect 
of  the  mean  service  time  (by  looking  at  the  figure's  columns),  the  effects  of  traffic  in- 
tensity (by  comparing  the  rows  in  a  sub-figure),  and  the  effect  of  the  combinations  of 
interarrival  and  service  time  distributions  (by  looking  at  the  columns  in  a  sub-figure). 
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However,  in  order  to  observe  the  effects  of  changes  to  the  probability  that  the  first  unit 
in  the  queue  selects  routel,  P,,  we  need  Figures  16  through  18  side-by-side. 

It  is  immediately  clear  in  Figures  16  through  18  that  traffic  intensity  p  has  a 
major  effect  on  the  mean  waiting  time  and,  in  fact,  it  is  known  from  theor}'  that  the 
mean  waiting  time  in  many  situations  is  proportional  to  the  reciprocal  of  (1-p),  [Ref  6: 
p.VIII.8].  Beyond  this,  it  appears  to  be  generally  true  that  the  regular-regular  combi- 
nation of  interarrival  and  service  times  at  the  right  most  position  in  the  six  skyscraper 
plots  of  Figure  16  is  best  for  traffic  performance;  on  the  other  hand,  the  skewed-skewed 
combination  at  the  left  most  position  in  the  six  plots  is  the  worst  case  for  all  values  of 
traffic  intensity.  In  addition,  as  the  mean  service  time  doubles  the  mean  waiting  time 
also  doubles.  Finally,  the  effects  of  increasing  mean  service  time  from  2.5  to  5.0  together 
with  other  factor  combinations  turns  out  to  be  consistent  with  different  values  of  P,. 

To  analyze  the  effects  of  the  factor  P,  requires  comparing  Figure  16  (Pi  =  0.75) 
through  Figure  18  (P,  =  0.25).  For  instance,  the  mean  waiting  time  of  the  east  direction 
is  almost  five  times  that  of  the  other  directions  in  Figure  16  with  P,  =  0.75.  On  the  other 
hand.  Figure  17  with  Pi  =  0.5  shows  that  the  mean  waiting  times  for  the  south  direction 
are  relatively  higher  than  the  north  and  east  directions.  This  is  because  the  frequency 
of  movements  on  each  route  varies  as  the  value  of  Pj  changes  in  conjunction  with  dif- 
ferent route  movement  conditions. 

Using  Figure  22  for  P,  =  1  and  traffic  intensity  =  0.5,  we  found  that  the  system 
is  stable  in  the  north  and  south  directions  but  can  be  unstable  in  the  east  direction  if 

mean  interarrival  time  of  east  units  .  ^ 


mean  interarrival  time  of  north  or  south  units 


A  more  sophisticated  analysis  of  the  results,  for  instance  the  test  to  see  if  there 
still  remains  a  difference  in  mean  waiting  times  of  three  directions  due  to  P,  after  ad- 
justing other  factors,  is  not  done  in  this  paper,  [Ref  6:  p. VIII. 9]. 
2.      Proposed  Road  System  (Model  2) 

The  results  of  the  simulation  for  Model  2  are  shown  in  Figures  19  through  21 
and  the  waiting  times  in  these  figures  are  compared  with  the  corresponding  values  of 
Figures  16  through  18.  The  effects  of  the  factors  found  in  the  first  model  are  consistent 
in  this  model.  We  can  expect  that  the  overall  traffic  performance  in  Model  2  would  in- 
crease because  the  new  system  allows  more  combinations  of  movements,  that  is,  the  re- 
strictions to  movements  are  reduced.  The  benefits  of  the  new  system  for  the  north 
direction  removes  two  possible  blocking  routes,  east  routel  and  south  routel,  which  are 
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blocking  units  in  the  first  model  for  north  route2.  As  for  the  south  direction,  east  route  1 
and  north  route2  no  longer  block  south  routel.  Finally  in  the  case  of  the  east  direction, 
east  routel  is  no  longer  blocked  by  south  routel  or  north  route2.  Thus,  overall  move- 
ment conditions  are  improved  and  waiting  times,  as  would  be  expected,  are  decreased  for 
the  three  directions  in  this  model. 

B.      COMPARISON  OF  THE  TWO  MODELS 

In  order  to  compare  the  results  of  the  two  simulations,  the  sample  means  and 
standard  deviations  of  the  first  500  average  (running)  waiting  times  for  each  direction 
(see  equations  3.7  and  3.9)  are  tabulated  in  Tables  19  through  21  of  Appendix  A,  so  that 
a  direct  comparison  of  the  waiting  times  for  the  two  models  is  possible  on  a  pair-wise 
basis.  The  fact  that  the  new  model  is  more  efficient  than  the  current  model  in  traffic 
performance  is  evident.  However,  the  criteria  to  prefer  the  new  system  to  the  old  could 
be  more  than  just  the  relative  differences  in  time.  Numbers  of  additional  units  deployed 
per  unit  of  time,  the  decrease  in  forces  susceptible  to  enemy  air  and  ground  ambushes, 
are  a  few  of  the  criteria  that  could  be  used  to  compare  the  two  systems. 
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Model  1   and  P.,  =  0.75  for  all  directions 
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Figure  16.  Effects  of  major  factors  for  Model  1:  Each  bar  in  the  cell  shows  the 
mean  estimate  of  the  first  500  (Nl  =  500)  waiting  times  for  each  di- 
rection's queue.  The  notation  SS,  SR,  RS,  and  RR  on  the  X  axis  rep- 
resent combinations  of  interarrival  and  service  time  distributions  where 
S  and  R  means  'skewed'  and  'regular'  interarrival  or  service  times  re- 
spectively. For  instance,  the  notation  SR  indicates  'skewed'  interarrival 
and  'regular'  service  time  distributions  were  used  for  the  simulation. 
These  are  achieved  by  using  different  Gamma  shape  parameters  A,,  /:2 
for  the  interarrival  and  service  times.  The  effects  of  the  traffic  intensity 
and  the  combination  of  interarrival  and  service  time  distributions  are 
obvious  and  consistent  in  the  six  skyscraper  plots.  The  mean  waiting 
times  in  the  east  direction  are  very  large  and  unstable  while  the  other 
two  directions  are  stable.  This  is  discussed  further  in  Figure  22.  Note 
that  the  scales  are  different  between  the  six  sub-figures. 
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Model  1  and  P^  =  0.50  for  all  directions 

Traffic  Intensity  =  0.5 

Kflean  Service  Time  =  2.5 
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Figure  17. 


Effects  of  probability  of  taking  route  1(0.50)  for  Model  I:  See  Figure 
16  for  the  explanation  of  SS,  SR,  RS,  and  RR  in  the  plot.  The  mean 
waiting  times  in  this  plot  correspond  to  a  mean  service  time  2.5  and 
traffic  intensity  of  0.5  for  each  direction.  The  effects  of  decreasing  P, 
from  0.75  in  Figure  16  to  0.50  in  this  figure  is  to  decrease  the  east's 
waiting  time  and  to  increase  the  south's  waiting  time.  The  decrease  in 
the  east's  waiting  time  is  because  the  probability  that  east  units  take 
east  route2  increases  to  0.50  from  0.25  in  Figure  16  and  it  is  easier  for 
east  units  to  move  to  east  route2  than  east  route  1  since  east  route  1  has 
three  blocking  routes.  On  the  other  hand,  there  is  only  one  blocking 
route  for  east  route2.  The  reason  for  the  increase  in  south's  waiting 
time  is  because  more  north  units  select  north  route2  and  more  east 
units  select  east  route2  which  combine  to  increase  the  difficulty  for 
south  units  to  move  under  Pi  =  0.50  vs  Pi  =0.75. 
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Model  1   and  P^  =  0.25  for  all  directions 

Traffic  Intensity  =  0.5 

Mean  Service  Time  =  2.5 
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Figure   18.      Effects  of  probability  of  taking  route  1(0.25)  for  Model  1:     See  Figure 

16  for  the  explanation  of  SS,  SR,  RS,  and  RR  in  the  plot.  The  mean 
waiting  times  in  this  figure  are  compared  with  the  corresponding  values 
in  the  third  row  of  the  upper  three  plots  of  Figure  16  and  against  Figure 
17.  The  major  effects  of  decreasing  the  value  of  Pj  from  0.50  in  Figure 

17  to  0.25  in  this  figure  is  to  decrease  east's  and  south's  waiting  times. 
The  decrease  in  east's  mean  waiting  time  is  because  now  east  units  se- 
lect east  route  2  with  probability  0.75  and  east  route2  is  only  blocked 
by  south  route  1,  which  is  selected  by  south  units  with  probability  0.25. 
That  is,  the  probability  that  an  arriving  east  unit  is  blocked  becomes 
lower  than  in  Figures  16  or  17.  The  decrease  in  the  south  waiting  times 
can  be  contributed  to  the  difference  in  the  probability  of  taking  each 
route  and  the  interactions  of  movements  on  the  routes. 
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Model  2  and  Pi  =  0.75  for  all  directions 
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Figure   19. 


Effects  of  major  factors  for  Model  2:  See  Figure  16  for  the  explanation 
of  SS,  SR,  RS,  and  RR  in  the  plots.  The  waiting  times  on  the  north, 
south  and  east  directions  decreased  almost  2,  6  and  more  than  10  times 
respectively  from  the  corresponding  waiting  times  of  Figure  16  in  the 
first  model.  The  difTerences  in  the  rates  are  because  the  performance 
benefits  of  the  new  system  are  different  for  each  direction.  For  example 
with  Pi  =  0.75,  the  benefits  of  new  east  routel  and  north  route2,  are 
different  since  approximately  75%  of  the  east  units  use  east  routel  and 
only  about  25%  of  the  north  units  use  north  route2.  However,  the 
overall  effects  of  traffic  intensity  and  mean  service  time  remain  the 
same  as  in  Model  1. 
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Model  2  and  P^  =  0.50  for  all  directions 

Traffic  Intensity  =  0.5 

Mean  Service  Time  =  2.5 
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Figure  20.  Effects  of  probability  of  taking  route  1(0.50)  for  Model  2:  See  Figure 
16  for  the  explanation  of  SS,  SR,  RS,  and  RR  in  the  plot.  The  mean 
waiting  times  in  this  figure  can  be  compared  with  the  corresponding 
values  in  Figure  17  to  see  Model  2's  effect.  The  corresponding  values 
in  the  third  row  of  each  direction  for  the  upper  three  plots  in  Figure  19 
depict  the  effects  of  changes  to  P,  in  Model  2.  Since  each  direction  has 
only  one  blocking  route  and  P,  =  0.5,  it  is  reasonable  to  expect  that  the 
mean  waiting  times  of  the  three  directions  are  approximately  equal. 
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Model  2  and  P^  =  0.25  for  all  directions 

Traffic  Intensity  =  0.5 

Mean  Service  Time  =  2.5 


ss 


(nterorrWoUSerV-ce 


RS 

-nme  Distribution 


Figure  21. 


Effects  of  probability  of  taking  route  1(0.25)  for  Model  2:  See  Figure 
16  for  the  explanation  of  SS,  SR,  RS,  and  RR  in  the  plot.  It  is  again 
obvious  that  the  mean  waiting  times  decrease  in  comparison  with  the 
corresponding  case  of  the  first  model  (Figure  18).  The  south  and  north 
have  longer  waiting  times  than  the  east  because  north  and  south  units 
select  north  route2  and  south  route2  with  probability  0.75.  On  the 
other  hand,  east  units  interact  only  with  the  north  and  south  units 
which  take  routes  blocking  east  units  with  probability  0.25. 
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Model  1   and  P^  =  0.75  for  all  directions 
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Figure  22.  Stability  changes  on  east  queue  (Model  1):  This  figure  shows  the 
specific  characteristics  of  the  intersection  given  iModel  1  with  P,  =  0.75 
and  traffic  intensity  =  0.5.  For  these  given  input  parameters,  this  figure 
shows  changes  of  stability  on  the  east  queue  given  changes  to  the  mean 
interarrival  time  ratios  for  the  three  directions.  The  changes  are  con- 
sistent with  any  combination  of  interarrival  and  service  time  distrib- 
utions, i.e.  SS,  SR,  RS,  and  RR  (see  Figure  16  for  an  explanation  of 
SS,  SR,  RS,  and  RR).  The  results  in  this  figure  are  from  skewed 
interarrival  and  service  time  distributions  (see  page  21).  The  numbers 
on  the  X  axis;  1.0,  1.1...  represent  the  ratio  of  (mean  interarrival  time 
of  east  units  /  mean  interarrival  time  of  north  or  south  units).  As  seen 
in  Figure  16,  the  east  queue  is  unstable  if  the  arrival  rate  from  the  east 
is  the  same  as  that  of  north  and  south,  and  becomes  stable  as  the  ar- 
rival rate  from  the  east  decreases  with  respect  to  the  north  or  south 
arrival  rate.  North's  and  south's  stability  was  not  afiected  by  the  ratio 
changes. 
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VI.      CONCLUSIONS  AND  RECOMMENDATIONS 

This  research  has  demonstrated  the  effects  of  traffic  conditions  on  traffic  flow  at  two 
kinds  of  three  way  intersections.  We  have  seen  how  traffic  along  the  roads  of  interest 
varies  as  the  conditions  change.  The  question  of  how  we  interpret  mathematical  terms 
such  as  Gamma  distributions  and  other  random  variables  in  the  real  world  still  remains 
uncertain.  That  is,  how  do  we  know  actual  interarrival  and  service  times  would  be  dis- 
tributed with  certain  characteristics?  For  example,  we  have  assumed  a  Gamma  distrib- 
ution for  interarrival  and  service  times  in  this  model,  however  this  assumption  may  be 
completely  different  from  the  real  world  situation. 

As  for  service  time,  it  seems  as  if  many  factors  affect  service  in  wartime  in  addition 
to  the  physical  restrictions  of  a  road  system.  As  mentioned  before,  major  problems  such 
as  refugee  flow  along  the  roads  should  be  studied  separately.  But  we  can  assume  a  cer- 
tain distribution  for  service  time  for  the  purpose  of  this  study. 

For  interarrival  times,  however,  we  can  collect  the  appropriate  data  by  conducting 
realistic  war  game  simulation.  As  the  war  game  situation  changes,  the  target  level  forces 
are  assumed  to  complete  their  movement  according  to  the  combat  model's  movement 
schedule,  [Ref  12:  p.  3.15-23].  The  movement  speed  is  assumed  to  vary  according  to  the 
characteristics  of  the  theater.  However,  the  current  models  do  not  consider  the  road 
movement  interaction  of  forces  which  may  create  critical  delays  in  combat  force  oper- 
ations. Therefore  we  must  consider  the  road  movement  interaction  of  combat  forces 
given  various  wartime  conditions  and  find  out  whether  the  movements  of  the  target  level 
forces  in  the  wargame  simulation  can  be  completed  in  a  timely  manner  in  the  real  world. 
This  study  enables  us  to  analyze  the  delays  in  combat  force  operations  by  considering 
the  interaction  of  convoys  which  have  to  traverse  an  intersection  of  interest.  If  we  create 
a  road  march  plan  for  the  forces  according  to  a  war  game  simulation  situation  then  we 
can  collect  the  expected  arrival  times  of  each  unit  at  the  intersection  of  interest.  The 
expected  arrival  times  at  an  intersection  of  interest  obtained  from  this  procedure  could 
be  developed  into  a  distribution  for  interarrival  times. 

Thus  we  can  apply  this  model  to  the  real  world  with  the  given  assumptions  and  can 
determine  expected  traffic  flows  in  the  area  of  interest.  In  this  thesis,  we  study  only  one 
three  way  intersection  in  an  attempt  to  describe  its  overall  traffic  performance.  It  would 
be  especially  useful  to  apply  this  model  to  a  three  way  intersection  which  is  believed  to 
be  a  traffic  choke  point  in  a  specific  area  of  interest.    Further  studies  about  more  ag- 
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gregated  road  systems  along  with  war  game  simulation  are  needed  to  further  understand 
the  potential  for  wartime  traffic  problems  in  the  Korean  theater. 
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APPENDIX  A.      SAMPLE  OUTPUT  OF  THE  SIMULATION  MODELS. 


Table  10.     SAMPLE   SECTION   OF  DIAGNOSTIC   PRINTOUT  OF  MODEL 

1:  Simulation  conditions;  Pj  =0.75,  mean  service  time  =2.5,  TI(trafric 
intensity)  =  0.5,  interarrival  and  service  time  distribution  =  SS  (see  Table 
8  and  Figure  16). 
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13.62 

.82 

1 

0 

2 

0 

3 

0 

0 

NNU 

14.58 

.  33 

1 

0 

3 

0 

3 

0 

0 

NS 

15.33 

.  26 

1 

0 

2 

0 

3 

0 

0 

NNU 

15.62 

.  08 

1 

0 

3 

0 

3 

0 

0 

NEU 

16.  04 

.  73 

1 

0 

3 

0 

4 

0 

0 

SS 

16.25 

.50 

0 

0 

3 

0 

4 

0 

0 

NNU 

16.32 

.20 

0 

0 

4 

0 

4 

0 

0 

NEU 

16.66 

.57 

0 

0 

4 

0 

5 

0 

0 

NS 

17.40 

.44 

0 

1 

4 

0 

4 

1 

0 

0 

ES 

17.92 

.91 

0 

2 

3 

0 

3 

0 

0 

0 

NS 

18.  13 

.31 

0 

2 

2 

0 

3 

0 

0 

0 

ES 

19.23 

.55 

0 

0 

2 

0 

3 

0 

0 

1 

NS 

19.96 

.91 

0 

0 

1 

2 

0 

2 

2 

0 

0 
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Table  11. 

SAMPLE  SECTION  OF  SINGLE.RUN.REPORT  OF  MODEL  1. 

AT 

arrival   time  of   a  unit. 

Ql 

queue   length  after   an  arrival. 

Q2 

queue   length  before  departure 

of  a  unit. 

WT 

waiting  time  of  each  unit. 

ST 

departure  time  of  a  unit. 

ET 

time  of  service  completion. 

DIR;    route  number  used  by  a  unit. 

NO 

AT              Ql            WT                   ST 

ET 

Q2 

DIR 

1 

.39            0            0.                       .39 

4.62 

0 

2 

.51             1            4.10              4.62 

8.  83 

1 

3 

7.51             1             1.32               8.83 

15.33 

1 

4 

9.61            1            5.72             15.33 

17.40 

3 

0 

r 

t 
h 

5 

13,62            2            4.30             17.92 

18.  13 

4 

6 

14.58            3            3.55             18.13 

19.96 

3 

7 

15.62            3            5.93            21.55 

22.  74 

3 

8 

16.32            4            6.42            22.74 

24.09 

2 

9 

20.48            3            5.54            26.02 

27.62 

2 

10 

24.35             2             3.26             27.62 

29.84 

1 

1 

8.62            0            0.                    8.62 

8.  96 

0 

2 

13.  17             0            0.                  13.  17 

16.  25 

0 

s 

3 

20.13             1            5.89            26.02 

28.  28 

2 

0 

1  -1 

4 

21.65            2            6.62            28.28 

28.54 

1 

5 

32.  76            0            0.                  32.  76 

38.  27 

0 

u 

6 

35.34             1            2.93            38.27 

38.66 

1 

t 
h 

7 

38.44             1            3.34            41.78 

42.83 

2 

8 

39.83            2            3.90            43.72 

44.42 

1 

9 

47.14             1               .32             47.47 

52.  37 

2 

10 

47.15            2            5.22            52.37 

52.  92 

3 

1 

7.08             1           10.32             17.40 

17.92 

5 

2 

7.88             2           10.04             17.92 

19.23 

4 

3 

8.79            3          11.16             19.96 

21.55 

3 

E 

4 

16.04            4            5.50            21.55 

22.85 

3 

a 

5 

16.66            5            6.19            22.85 

23.57 

2 

s 

6 

20.  12            3            3.  96            24.  09 

26.02 

2 

t 

7 

22.94            2            6.90            29.84 

30.  70 

3 

8 

28.43            2            9.84            38.27 

38.54 

5 

9 

28.69            3            9.85            38.54 

41.  78 

4 

10 

30.91             3           11.92             42.83 

43.  72 

5 
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Table  12. 

SAMPLE  SECTION  OF  SUM.REPORT  OF  MODEL  1 

• 

AQn:    average 

queue 

length( 

running)    at 

each  north 

unit' s 

departure. 

AWn:    average 

waiting  time( 

running)    at 

each  north 

unit's 

departure. 

AQs   and 

AWs 

are   for 

south 

units ,   AQe 

and  AWe  are 

for  east  units.              | 

No 

AQn 

AWn 

AQs 

AWs 

AQe 

AWe 

1 

0. 

0. 

0. 

0. 

1.76 

10.32 

2 

.89 

2.05 

0. 

0. 

1.82 

10.  18 

3 

.61 

1.81 

.39 

1.96 

1.94 

10.51 

4 

.89 

2.  79 

.44 

3.  13 

2.01 

9.26 

5 

1.27 

3.09 

.38 

2.50 

2.01 

8.64 

6 

1.29 

3.  17 

.40 

2.57 

2.01 

7.86 

7 

1.45 

3.56 

.50 

2.68 

1.90 

7.  73 

8 

1.48 

3.92 

.52 

2.83 

2.32 

7.99 

9 

1.48 

4.  10 

.49 

2.56 

2.  33 

8.20 

10 

1.45 

4.01 

.65 

2.82 

2.47 

8.57 

11 

1.  34 

3.74 

.66 

3.04 

2.53 

8.75 

12 

.86 

3.42 

.67 

2.97 

3.05 

9.95 

13 

.  79 

3.  16 

.85 

3.20 

3.47 

11.42 

14 

.  78 

2.97 

.91 

3.43 

3.60 

12.85 

15 

.  71 

2.90 

1.23 

4.  10 

3.93 

14.53 

16 

.72 

2.93 

1.43 

5.04 

3.94 

15.93 

17 

.67 

2.88 

1.44 

5.60 

3.95 

16.96 

18 

.62 

2.  72 

1.52 

6.36 

3.99 

17.54 

19 

.69 

2.92 

1.52 

6.66 

4.  18 

17.84 

20 

.  73 

3.09 

1.55 

6.69 

4.20 

18.04 

21 

.  74 

3.  19 

1.63 

6.88 

4.21 

17.77 

22 

.  72 

3.04 

1.64 

6.96 

4.51 

17.83 

23 

.63 

2.91 

1.61 

7.36 

4.56 

17.92 

24 

.58 

2.  79 

1.60 

7.  16 

4.63 

18.20 

25 

.57 

2.68 

1.61 

7.04 

5.03 

18.  73 

26 

.55 

2.65 

1.60 

6.94 

5.  11 

19.  13 

27 

.57 

2.76 

1.31 

6.69 

5.  18 

19.54 

28 

.61 

2.84 

1.27 

6.47 

5.40 

20.  12 

29 

.63 

3.05 

1.24 

6.45 

5.61 

20.87 

30 

.63 

3.03 

1.54 

6.96 

5.72 

21.29 
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Table  13.     SAMPLE  SECTION   OF   DIAGNOSTIC   PRINTOUT  OF  MODEL 

2.:  Simulation  conditions;  Pi  =0.75,  mean  service  time  =2.5,  Tl(trafric 
intensity)  =  0.5,  interarrival  and  service  time  distribution  =  SS  (see  Table 
8  and  Figure  16). 

NNU:  new. north. unit ,  NSU:  new. south,  unit ,  NEU:  new.  east,  unit 

NS:  north. service,    SS:  south. service,    ES:  east. service 

U:  a  uniform  random  variable  used  for  the  next  route  selection. 

Movements:  movements  at  the  intersection  after  the  left  event. 

N,S,E:  north,  south,  and  east  direction. 

Queue:  number  of  units  in  each  queue  after  the  left  event. 

Next  route:  routes  selected  for  the  first  unit  in  each  queue. 


Event 


Time 


Movements 

Queue 

Next  route 

u 

(N 

,  s 

,  E) 

(N,S 

,E) 

(N,S 

,E) 

45 

0 

0 

0 

0 

0 

0 

0 

0 

43 

0 

0 

1 

0 

0 

0 

0 

0 

66 

0 

0 

0 

0 

0 

0 

0 

0 

52 

0 

0 

0 

0 

1 

0 

0 

04 

0 

0 

0 

1 

0 

0 

44 

0 

0 

0 

2 

0 

0 

34 

1 

0 

0 

2 

0 

0 

54 

1 

0 

0 

3 

0 

0 

42 

0 

1 

0 

2 

0 

0 

23 

0 

0 

0 

2 

0 

0 

54 

0 

0 

2 

0 

2 

0 

0 

07 

0 

1 

2 

0 

2 

0 

0 

82 

0 

1 

3 

0 

2 

0 

0 

33 

0 

1 

4 

0 

2 

0 

0 

26 

1 

0 

3 

0 

2 

0 

0 

1 

08 

1 

0 

4 

0 

2 

0 

0 

1 

73 

1 

0 

4 

0 

3 

0 

0 

1 

50 

0 

0 

4 

0 

3 

0 

0 

1 

20 

0 

0 

5 

0 

3 

0 

0 

1 

57 

0 

0 

5 

0 

4 

0 

0 

1 

44 

0 

1 

5 

0 

3 

1 

0 

0 

91 

0 

2 

4 

0 

2 

0 

0 

0 

31 

0 

2 

3 

0 

2 

0 

0 

0 

55 

0 

0 

3 

0 

2 

0 

0 

1 

91 

2 

0 

2 

0 

1 

0 

0 

0 

05 

2 

0 

2 

0 

2 

0 

0 

0 

37 

2 

1 

2 

0 

2 

0 

0 

0 

64 

2 

1 

3 

0 

2 

0 

0 

0 

99 

2 

1 

2 

0 

2 

0 

0 

0 

55 

2 

0 

2 

0 

2 

0 

0 

0 

84 

2 

0 

2 

2 

0 

1 

0 

0 

0 

83 

2 

0 

2 

2 

1 

1 

0 

2 

0 

NNU 

.  39 

NNU 

.51 

NS 

4.62 

NEU 

7.08 

NNU 

7.51 

NEU 

7.88 

NSU 

8.  62 

NEU 

8.  79 

NS 

8.83 

SS 

8.96 

NNU 

9.61 

NSU 

13.  17 

NNU 

13.62 

NNU 

14.58 

ES 

15.33 

NNU 

15.62 

NEU 

16.04 

SS 

16.  25 

NNU 

16.  32 

NEU 

16.66 

NS 

17.40 

ES 

17.92 

NS 

18.  13 

ES 

19.23 

NS 

19.96 

NEU 

20.  12 

NSU 

20.  13 

NNU 

20.48 

NS 

21.  14 

SS 

21.43 

ES 

21.55 

NSU 

21.65 
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Table  14. 

SAMPLE  SECTION  OF  SINGLE.RUN.REPORT  OF  MODEL  2. 

AT:      arrival   time  of  a  unit. 

Ql 

queue   length  after   an  arrival. 

Q2 

queue   length  before  departure 

of  a  unit. 

WT 

waiting  time  of   each  unit. 

ST 

departure  time  of  a  unit. 

ET 

time  of  service  completion. 

DIR:    route  number  used  by  a  unit. 

NO 

AT              Ql            WT                  ST 

ET 

Q2 

DIR 

1 

.39            0            0.                       .39 

4.62 

0 

2 

.51             1            4.10               4.62 

8.83 

1 

N 

3 

7.51             1             7.82             15.33 

17.40 

4 

4 

9.61            2            8.31             17.92 

18.  13 

5 

0 

5 

13.62            3            4.51             18.13 

19.96 

4 

t 
h 

6 

14.58            4            5.37             19.96 

21.  14 

3 

7 

15.62            4            5.52            21.14 

22.49 

3 

8 

16.32            5             7.89             24.21 

26.46 

2 

9 

20.48             3            8.21            28.69 

28.95 

2 

10 

24.35             2            5.47            29.82 

34.99 

2 

1 

8.62            0            0.                    8.62 

8.96 

0 

2 

13.  17             0            0.                  13.  17 

16.25 

0 

s 

3 

20.  13            0            0.                  20.  13 

21.43 

0 

0 

4 

21.65             1               .84            22.49 

24.09 

1 

5 

32.  76            0            0.                  32.  76 

38.27 

0 

u 

6 

35.34             1             2.93             38.27 

39.  32 

1 

t 

1 

7 

38.44             1               .88             39.32 

40.01 

1 

h 

8 

39.83             1                .19             40.01 

44.92 

1 

9 

47.  14            0            0.                  47.  14 

47.69 

0 

10 

47.15             1               .54            47.69 

49.  75 

2 

1 

7.08             1             1.75               8.83 

15.33 

3 

2 

7.88            2            9.52             17.40 

17.92 

4 

3 

8.79             3            9.13             17.92 

19.23 

3 

E 

4 

16.04            3            3.91             19.96 

21.55 

2 

a 

5 

16.66            4            4.89             21.55 

22.27 

2 

s 

6 

20.12            2            2.15             22.27 

24.21 

1 

t 

7 

22.94            1            3.52            26.46 

28.69 

1 

8 

28.43             1               .53            28.95 

29.82 

2 

9 

28.69             2            6.29             34.99 

35.26 

4 

10 

30.91            2            4.75             35.66 

38.90 

3 

60 


Table  15.     SAMPLE  SECTION  OF  SUM.REPORT  OF  MODEL  2. 


AQn:  average 

queue 

length( 

running)  at 

each  north 

unit' s 

departure. 

AWn:  average 

waiting  time( 

running)  at 

each  north 

unit ' s 

departure. 

AQs  and  AWs 

are  for 

south 

units,  AQe 

and  AWe  are 

for  east  units.      | 

No 

AQn 

AWn 

AQs 

AWs 

AQe 

AWe 

1 

0. 

0. 

0. 

0. 

.31 

1.  75 

2 

.89 

2.05 

0. 

0. 

1.26 

5.63 

3 

1.31 

3.97 

0. 

0. 

1.31 

6.80 

4 

1.77 

5.06 

.04 

.21 

1.  38 

6.08 

5 

1.80 

4.95 

.03 

.  17 

1.42 

5.  84 

6 

1.91 

5.02 

.  10 

.63 

1.41 

5.22 

7 

1.95 

5.09 

.  12 

.66 

1.32 

4.98 

8 

1.95 

5.44 

.  12 

.60 

1.23 

4.42 

9 

1.95 

5.  75 

.10 

.54 

1.34 

4.63 

10 

1.92 

5.  72 

.  11 

.54 

1.38 

4.64 

11 

1.  78 

5.70 

.  15 

.68 

1.43 

4.68 

12 

1.31 

5.23 

.21 

.97 

1.42 

4.  70 

13 

1.22 

5.00 

.22 

.99 

1.39 

4.56 

14 

1.21 

5.02 

.35 

1.  18 

1.39 

4.31 

15 

1.  15 

4.  77 

.54 

1.63 

1.28 

4.02 

16 

1.  15 

4.58 

.  73 

2.  37 

1.22 

4.00 

17 

1.03 

4.31 

.  74 

2.  73 

1.  18 

3.91 

18 

.92 

4.07 

.  79 

3.  19 

1.  08 

3.  73 

19 

.91 

3.86 

.79 

3.  23 

.  93 

3.53 

20 

.88 

3.  78 

.74 

3.06 

.90 

3.39 

21 

.88 

3.67 

.  72 

2.94 

.80 

3.23 

22 

.83 

3.69 

.69 

2.  80 

.78 

3.08 

23 

.77 

3.53 

.68 

2.68 

.78 

2.99 

24 

.72 

3.53 

.59 

2.57 

.90 

3.25 

25 

.73 

3.55 

.58 

2.47 

.92 

3.41 

26 

.74 

3.55 

.58 

2.51 

1.01 

3.53 

27 

.74 

3.52 

.47 

2.41 

1.  04 

3.66 

28 

.77 

3.61 

.47 

2.39 

1.  11 

3.97 

29 

.77 

3.  73 

.46 

2.31 

1.  11 

4.29 

30 

.98 

4.  23 

.43 

2.23 

1.  12 

4.22 
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Table  16.     SIMULATION    OUTPUTS     FOR    VARIOUS    CONDITIONS    OF 
MODEL  1 


P^:      probability  taking 

routel. 

MS:      mean  service  time. 

TI:      traffic   intensity(see  Table  8). 

A. S:    interarrival   and  service  time  distribution  (see  Figure   16). 

WT:      mean  waiting  time   ( 

;see  equation  (3. 7)    in  chapter   III. ) 

SD:       std.  dev.    of  WT(see 

equation   (3.8)    in  chapter   III.) 

Pi 

MS 

TI 

A.S 

NORTH                     SOUTH                     EAST 
WT            SD            WT            SD            WT            SD 

SS 

16.  71 

1.63      28.40        2.  39    155.  93      18.  38 

0.5 

SR 

11.38 

.  77      22.67        1.82    120.57      10.  39 

RS 

10.  90 

.65      26.29        3.21    159.92      16.34 

RR 

4.  87 

.29      13.31           .96      99.75      12.45 

SS 

7.27 

.39        9.72          .75      16.41        1.35 

2.5 

0.4 

SR 
RS 

4.93 
4.  72 

.32        6.23          .37      11.78        1.36 
.33        6.90          .38      11.50          .66 

RR 

2.09 

.08        3.18          .15        4.53          .18 

SS 

4.49 

.20        5.66          .22        8.54          .67 

0.33 

SR 

2.68 

.10        3.30          .10        4.41          .19 

RS 

2.45 

.12        2.98          .10        4.09          .18 

RR 

1.  15 

.04        1.49          .05        1.92          .04 

0.75 

SS 

37.56 

2.83      82.31      13.  74   337.  10      27.53 

0.5 

SR 

23.54 

1.47      50.80        5.42   258.47      24.42 

RS 

22.59 

1.94     50.  35        5.91   344.53      37.02 

RR 

10.40 

.49      30.51        2.63    163.22      20.  72 

SS 

14.46 

.96      16.68          .81     32.62        3.86 

5.0 

0.4 

SR 
RS 

9.  11 
9.00 

.38      12.22          .53      22.31        1.47 
.46      12.77        1.07      22.56        1.79 

RR 

3.90 

.11        5.91          .27        7.99          .30 

SS 

9.84 

.49      10.49          .43      16.71        1.17 

0.33 

SR 

5.58 

.17        6.27          .24        8.77          .36 

RS 

5.39 

.25        6.54          .35        8.86          .46 

RR 

2.43 

.08        3.00          .14        3.83          .15 
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Table  17. 

CONTINUED  FROM  THE  PREVIOUS  TABLE. 

Pi 

MS 

TI 

A.S 

NORTH 
WT            SD 

SOUTH 
WT            SD 

EAST 
WT            SD 

SS 

21.13        2.47 

67.87      12.09 

17.  75 

1.  30 

0.5 

SR 
RS 
RR 

11.42        1.06 
14.80        1.57 

29.23        3.71 
47.61        8.42 

9.33 
11.88 

.66 

1.  64 

5.40          .33 

16.00        1.97 

4.50 

.20 

SS 

6.85          .30 

10.05        1.18 

7.  19 

.43 

2.5 

0.4 

SR 
RS 
RR 

4.49          .20 
4.51          .21 

6.41          .47 
7.31          .41 

4.43 
4.  14 

.21 
.25 

1.87          .06 

3.20          .18 

1.73 

.06 

SS 

4.37          .19 

4.99          .27 

4.  16 

.  15 

0.33 

SR 
RS 
RR 

2.59          .11 
2.64          .13 

3.04          .16 
3.34          .17 

2.57 
2.  38 

.  13 
.  13 

0.50 

1.11          .03 

1.27          .04 

.98 

.03 

SS 

36.04        3.69 

106.57      22.53 

39.94 

4.55 

0.5 

SR 
RS 
RR 

21.93        1.64 
25.03        2.11 

59.29        6.06 
79.50      12.27 

20.82 
26.06 

1.53 
1.86 

11.87          .89 

32.57        3.47 

9.93 

.64 

SS 

15.11        1.41 

21.31        1.41 

14.39 

.82 

5.0 

0.4 

SR 
RS 
RR 

8.30          .30 
9.03          .44 

12.25          .83 
14.88        1.01 

8.20 
8.50 

.57 
.39 

3.94          .12 

6.04          .29 

3.59 

.  14 

SS 

8.87          .62 

10.31          .47 

9.07 

.45 

0.33 

SR 
RS 
RR 

5.44          .21 
5.15          .23 

6.12          .34 
6.50          .25 

5.06 
4.61 

.  17 
.21 

2.24          .07 

2.91          .10 

2.06 

.06 
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Table  18. 

CONTINUED  FROM  THE  PREVIOUS  TABLE 

• 

Pi 

MS 

TI 

A.S 

NORTH 
WT            SD 

SOUTH 
WT            SD 

EAST 
WT 

SD 

SS 

16.77        1.21 

27.61        2.99 

6.  15 

.  39 

0.5 

SR 
RS 

12.16        1.34 
13.12          .98 

22.73        3.41 
28.20        5.13 

3.48 
3.44 

.  16 
.22 

RR 

5.79          .35 

11.00        1.36 

1.37 

.05 

SS 

7.44          .35 

8.53          .62 

3.29 

.  16 

2.5 

0.4 

SR 
RS 

4.58          .18 
4.  36          . 19 

4.  85          .  25 
5.20          .25 

2.06 
1.84 

.09 
.07 

RR 

1.95          .06 

2.12          .11 

.76 

.04 

SS 

4.53          .28 

4.61          .19 

2.43 

.09 

0.33 

SR 
RS 

2.53          .07 
2.49          .10 

2.71          .10 
2.87          .16 

1.47 
1.23 

.07 
.07 

0.25 

RR 

1.06          .03 

1.14          .05 

.48 

.02 

SS 

38.06        3.29 

76.42      11.18 

11.85 

.80 

0.5 

SR 
RS 

21.42        1.34 
25.02        2.25 

34.54        2.93 
50.95        6.00 

7.  16 
7.01 

.38 
.39 

RR 

10.98          .61 

17.03        1.16 

2.58 

.  11 

SS 

14.03          .82 

17.44        1.41 

6.97 

.40 

5.0 

0.4 

SR 
RS 

8.47          .28 
8.50          .51 

8.72          .37 
10.00          .47 

4.42 
4.04 

.30 
.23 

RR 

3.  90          .  16 

4.35          .21 

1.48 

.04 

SS 

8.91          .36 

8.26          .25 

4.  71 

.21 

0.33 

SR 
RS 

5.29          .17 
4.72          .13 

5.67          .21 
5.22          .19 

2.94 
2.  26 

.  10 
.08 

RR 

2.24          .08 

2.  34          .  06 

1.00 

.  03 
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Table  19. 

SIMULATION    OUTPUTS    FOR    VARIOUS 

CONDITIONS    OF 

MODEL  2 

Pj^:      probability  taking  routel. 

MS:      mean  service  time. 

TI:      traffic   intensity(see  Table  8). 

A. S:    interarrival   and  service  time  distribution  (see  Figure   16). 

WT:      mean  waiting  time   (see  equation  (3.7)    in 

chapter   III. ) 

SD:      std.  dev.    of  WT(see  equation  (3.8)    in  chapter   III.) 

Pi 

MS 

TI 

A.S 

NORTH                     SOUiH 
WT            SD            WT            SD 

EAST 
WT            SD 

SS 

9.09          .82        4.10          .21 

11.89        1.23 

0.5 

SR 

5.22          .27        2.61          .11 

6.66          .44 

RS 

6.14          .33        2.13          .08 

7.19          .46 

RR 

2.36          .13          .84          .03 

3.00          .18 

SS 

4.65          .26        2.75          .16 

5.32          .31 

2.5 

0.4 

SR 
RS 

2.87          .12        1.56          .05 
2.57          .10        1.26          .06 

2.96          .13 
2.95          .12 

RR 

1.07          .04          .51          .02 

1.19          .04 

SS 

3.03          .16        1.69          .07 

3.30          .14 

0.33 

SR 
RS 

1.87          .07        1.11          .04 
1.51          .06          .80          .04 

2.  11          .  12 
1.73          .09 

RR 

.63          .02          .27          .01 

.73          .02 

0.75 

SS 

19.81        1.66        8.38          .35 

26.59        1.92 

0.5 

SR 
RS 

11.62          .71        5.48          .48 
13.17          .81        4.74          .27 

15.15        1.44 
16.95        2.65 

RR 

4.75          .24        1.68          .06 

5.99          .32 

SS 

8.76          .48        5.18          .19 

10.18          .55 

5.0 

0.4 

SR 
RS 

5.67          .27        3.28          .09 
5.10          .23        2.31          .09 

6.39          .26 
5.85          .35 

RR 

2.08          .07          .84          .03 

2.40          .08 

SS 

5.71          .29        3.63          .17 

6.53          .24 

0.33 

SR 
RS 

3.91          .20        2.24          .07 
3.18          .13        1.52          .05 

4.15          .21 
3.44          .13 

RR 

1.27          .05          .57          .02 

1.50          .04 
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Table  20. 

CONTINUED  FROM  THE  PREVIOUS  TABLE 

• 

Pi 

MS 

TI 

A.S 

NORTH 
WT            SD 

SOUTH 
WT            SD 

EAST 
WT 

SD 

SS 

6.42          .46 

6.29          .34 

5.66 

.24 

0.5 

SR 
RS 

3.95          .19 
3.61          .19 

4.03          .20 
3.79          .20 

3.50 
3.96 

.09 
.  19 

RR 

1.40          .05 

1.44          .05 

1.39 

.04 

SS 

3.60          .19 

3.89          .19 

3.60 

.  14 

2.5 

0.4 

SR 
RS 

2.22          .07 
1.96          .09 

2.08          .07 
1.83          .08 

2.09 
1.95 

.08 
.  12 

RR 

.74          .02 

.74          .02 

.74 

.02 

SS 

2.35          .11 

2.28          .10 

2.40 

.09 

0.33 

SR 
RS 

1.49          .04 
1.15          .04 

1.49          .06 
1.18          .04 

1.43 
1.  18 

.05 
.07 

0.50 

RR 

.46          .01 

.49          .02 

.46 

.01 

SS 

13.03          .77 

14.24        1.06 

14.05 

1.  13 

0.5 

SR 
RS 
RR 

8.25          .48 
7.94          .49 

7.54          .26 
7.88          .47 

7.49 
8.61 

.36 
.51 

2.91          .13 

2.89          .13 

2.81 

.  12 

SS 

7.27          .28 

6.80          .30 

6.59 

.21 

5.0 

0.4 

SR 
RS 

4.  20          .  18 
3.39          .11 

4.45          .16 
3.45          .13 

4.58 
3.  73 

.  15 
.  13 

RR 

1.47          .05 

1.42          .03 

1.45 

.05 

SS 

4.58          .18 

4.43          .17 

5.  19 

.23 

0.33 

SR 
RS 

2.96          .12 
2.33          .08 

2.86          .07 
2.39          .10 

3.07 
2.45 

.  12 
.  11 

RR 

1.01          .04 

.96          .03 

.98 

.03 
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Table  21. 

CONTINUED  FROM  THE  PREVIOUS  TABLE 

• 

P. 

MS 

TI 

A.S 

NORTH 
WT            SD 

SOUTH 
WT            SD 

EAST 
WT 

SD 

SS 

11.13          .84 

12.27          .66 

4.51 

.24 

0.5 

SR 
RS 

6.38          .63 
5.64          .36 

7.24          .54 
6.69          .46 

2.59 
2.22 

.  10 
.  13 

RR 

2.31          .07 

2.84          .10 

.80 

.03 

SS 

4.  65          . 28 

4.81          .27 

2.24 

.  13 

2.5 

0.4 

SR 
RS 
RR 

2.65          .09 
2.62          .12 

2.81          .11 
3.01          .16 

1.67 
1.20 

.08 
.05 

1.01          .02 

1.12          .02 

.45 

.02 

SS 

2.85          .15 

3.14          .11 

1.  75 

.  09 

0.33 

SR 
RS 
RR 

1.87          .09 
1.60          .04 

1.91          .06 
1.85          .08 

1.  16 
.  76 

.03 

.04 

0.25 

.61          .02 

.75          .02 

.28 

.  01 

SS 

18.62        2.82 

20.37        1.70 

7.83 

.42 

0.5 

SR 
RS 
RR 

12.13        1.09 
12.60          .78 

13.01          .89 
13.08          .60 

5.27 
4.87 

.  30 
.24 

4.86          .26 

6.16          .26 

1.80 

.08 

SS 

8.54          .41 

9.02          .37 

5.  71 

.26 

5.0 

0.4 

SR 
RS 
RR 

5.76          .25 
4.98          .15 

6.13          .25 
5.24          .21 

3.  07 
2.37 

.  10 
.  12 

2.18          .05 

2.46          .08 

.91 

.  03 

SS 

6.  09          .  26 

6.29          .18 

3.64 

.  15 

0.33 

SR 
RS 
RR 

3.62          .13 
2.  99          .  08 

4.03          .15 
3.41          .10 

2.23 
1.48 

.  07 
.06 

1.27           .05 

1.47          .03 

.57 

.  02 
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APPENDIX  B.      A  PROGRAM  LISTING  OF  MODEL  1. 

PREAMBLE    ''an  intersection  simulation 
normally,  mode  is  undefined 
event  notices  include  north. unit,  generator, 

south. unit,  generator, 

east. unit,  generator, 

new. north,  unit , 
new. south,  unit, 
new.  east,  unit, 

north. service, 
south, service, 
east. service 

every  north. service  has  an  north. pointer 
every  south. service  has  an  south. pointer 
every  east. service  has  an  east. pointer 

define  north,  pointer, 
south. pointer, 
east. pointer  as  pointer  variables 

temporary  entities 

every  unit  may  belong  to  the  north. queue 

may  belong  to  the  south. queue 

may  belong  to  the  east. queue 

the  system  owns  the  north. queue  , 

the  south. queue, 
the  east. queue 

define  north,  queue, 
south,  queue 
and  east. queue  as  f ifo  sets 

define  north,  run, 
south,  run, 
east,  run, 
next. north,  route, 
next. south,  route, 
next. east,  route, 
Nl, 

replication,NR, 
in , j  n , kn , is , j  s , ks , ie , j  e , ke , 
saveseedl,saveseed2,saveseed3  as  integer  variables 

define  AQ.north.nl,  AQ.south.nl,  AQ.east.nl, 
AW.north.nl,  AW.south.nl,  AW.east.nl, 
alphal,  betal,  alpha2,  beta2, 
PI,  current,  time, 
sl,s2,s3,rl,r2,r3  as  real  variables 
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define  north. dir, south. dir,east. dir 

as  1-dimensional  integer  array 

define  arrival,  time.  of.  north,  unit, 

departure,  time.  of.  north,  unit, 
end.  time.  of.  north,  service, 
waiting,  time.  of.  north,  unit, 
north. queue. lengthl, 
north. queue. Iength2, 

arrival,  time.  of.  south,  unit, 
departure,  time.  of.  south,  unit, 
end.  time.  of.  south,  service, 
waiting,  time.  of.  south,  unit, 
south. queue,  lengthl, 
south,  queue.  Iength2, 

arrival,  time.  of.  east,  unit, 

departure,  time.  of.  east,  unit, 

end.  time.  of.  east,  service, 

waiting,  time.  of.  east,  unit, 

east,  queue. lengthl, 

and  east. queue. Iength2 , 

AQ.  of .  north,  unit ,  AQ.  of.  south,  unit ,  AQ.  of .  east,  unit , 

AW.  of .  north,  unit ,  AW.  of .  south,  unit ,  AW.  of .  east,  unit 

as  1-dimensional  real  array 

tally  ave.que.north.nl  as  the  mean  and 

sd.  ave.  que.  north,  nl  as  the  std.  dev  of  AQ.  north,  nl 
tally  ave.wt.north.nl  as  the  mean  and 

sd.ave.wt.north.nl  as  the  std.  dev  of  AW.north.nl 

tally  ave.que.south.nl  as  the  mean  and 

sd.  ave.  que.  south,  nl  as  the  std.  dev  of  AQ.  south,  nl 
tally  ave.wt.south.nl  as  the  mean  and 

sd.ave.wt.south.nl  as  the  std.  dev  of  AW.south.nl 

tally  ave. que. east. nl  as  the  mean  and 

sd.  ave.  que.  east,  nl  as  the  std.  dev  of  AQ.  east,  nl 
tally  ave.wt.east.nl  as  the  mean  and 

sd.ave.wt.east.nl  as  the  std.  dev  of  AW.east.nl 

accumulate  AQ. north. unit  as  the  mean  of  n. north. queue 
accumulate  AQ. south. unit  as  the  mean  of  n. south. queue 
accumulate  AQ. east. unit  as  the  mean  of  n. east. queue 

priority  order  is  north. unit. generator , 

south. unit. generator, 
east. unit,  generator 

END  ' ' preamble 


MAIN  ' '  program  starts  execution  here, 
define  s. que. north. 
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s.  wt.  north, 

s.  que.  south, 

s. wt. south, 

s.  que.  east, 

s.  wt.  east, 

inverse, traffic,  intensity  as  real  variables 


call  read. data 
saveseedl 
saveseed2 
saveseed3 


seed.  v(sl) 
seed. v(s2) 
seed.  v(s3) 


for  PI  =  0.25  to  0.75  by  0.25 

do 
also  for  alpha2  =  2.  5  to  5.  0  by  2.5 

do 
also  for  inverse  =  2  to  3  by  .  5 

do 

traffic,  intensity  =  1/inverse 

alphal  =  alpha2"inverse 

rl  =  1 

r2  =  1 

r3  =  1 
also  for  betal  =  1  to  4  by  3 

do 
also  for  beta2  =  1  to  4  by  3 

do 

print  2  line  thus 

PI    al    bl    a2    b2   TI 
print  2  line  with  PI , alphal, betal, alpha2, 

beta2, traffic. intensity  thus 


hJ«     W«^'^     W^^d     mjf  ^J^mJ^ 


for  replication  =  1  to  NR 
do 

let  seed. v(sl)  =  saveseedl 

let  seed. v(s2)  =  saveseed2 

let  seed. v(s3)  =  saveseed3 

call  initialize 
start  simulation 
call  single. run. report 
call  sum. report 


AQ.  north,  nl 
AQ.  south,  nl 
AQ.  east,  nl 
AW. north,  nl 
AW.  south,  nl 
AW.  east,  nl 


AQ.  of.  north.  unit(Nl) 
AQ.  of.  south.  unit(Nl) 
AQ.  of.  east.  unit(Nl) 
AW.  of.  north.  unit(Nl) 
AW.  of.  south.  unit(Nl) 
AW.  of.  east.unit(Nl) 


let  time,  v  =  0 

let  saveseedl  =  seed. v(sl) 

let  saveseed2  =  seed. v(s2) 
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let  saveseed3  =  seed. v(s3) 

let  in  =  0 
is  =  0 
ie  =  0 

jn  =  0 
js  =  0 
je  =  0 

kn  =  0 
ks  =  0 
ke  =  0 

if  replication  =  NR  and  ^fR  >  1 

s.  que.  north  =  sd.  ave.  que.  north,  nl/sqrt.  f(NR  -  1) 
s.  que.  south  =  sd.  ave.  que.  south,  nl/sqrt.  f(NR  -  1) 
s.  que.  east  =  sd.  ave.  que.  east,  nl/sqrt.  f(NR  -  1) 

s.wt.  north  =  sd.  ave.  wt.  north,  nl/sqrt.  f(NR  -  1) 
s.wt.  south  =  sd.  ave.  wt.  south,  nl/sqrt.  f(NR  -  1) 
s.wt.  east  =  sd.  ave.  wt.  east,  nl/sqrt.  f(NR  -  1) 


''print  estimators  of  mean  queue  length  (running) 
''at  the  Nlth  unit's  departure  and  mean  waiting  time 
''(running)  of  the  Nlth  unit  and  std. dev  of  each. 

print  5  lines  with  ave.que.north.nl, 

s.  que.  north, 
ave.  que.  south,  nl, 
s.  que.  south, 
ave.  que.  east,  nl, 
s. que. east, 

ave.  wt.  north,  nl  ,s.  wt.  north, 
ave.  wt.  south.  nl,s.  wt.  south, 
ave.  wt.  east,  nl  ,s.  wt.  east  thus 


See  main  routine  for  information  of  data. 


always 

release  arrival,  time.  of.  north,  unit , 

departure,  time.  of.  north,  unit , 
end.  time.  of.  north,  service, 
waiting,  time.  of.  north,  unit, 
north. queue. Iength2, 
arrival,  time.  of.  south,  unit, 
departure,  time.  of.  south,  unit, 
end.  time.  of.  south,  service, 
waiting,  t.-^me.  of.  south,  unit, 
south. queue. lengthZ, 
arrival,  time.  of.  east,  unit, 
departure,  time.  of.  east,  unit. 
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end.  time.  of.  east,  service, 

waiting,  time.  of.  east,  unit, 

east. queue.  Iength2, 

AQ.  of.  north,  unit , 

AQ.  of.  south,  unit, 

AQ.  of.  east,  unit, 

AW.  of.  north,  unit , 

AW.  of.  south,  unit , 

AW.  of.  east,  unit, 

north. queue. lengthl, 

south,  queue,  lengthl, 

east. queue,  lengthl, 

north. dir, south.  dir,east.  dir 


loop 


reset  the  totals  of  AQ. north. N1,AQ. south. N1,AQ.  east.  Nl, 

AW.  north.  Nl , AW.  south.  Nl , AW.  east.  Nl 
loop 
END  ^ 'MAIN 


ROUTINE  to  READ.  DATA 
print  4  lines  thus 


See  MAIN  routine  to  see  other  factors 
that  are  not  asked  here  ! 

print  2  lines  thus 

enter  the  seed  stream  numberl,2,3 
si,  s2  ,s3  (1  -  10) 

read  sl,s2,s3 

print  3  lines  thus 

enter  the  desired  number  of  units 
and  replications  simulated  ! 

Nl,  NR 
read  N1,NR 
END  ' '  READ.  DATA 


ROUTINE  to  INITIALIZE 

let  seed. v(sl)  =  saveseedl 
let  seed. v(s2)  =  saveseed2 
let  seed. v(s3)  =  saveseed3 

schedule  a  north. unit. generator  in 

gamma.  f(rl'<"alphal,betal,sl)  minutes 
schedule  a  south. unit. generator  in 

gamma. f(r2*alphal,betal,sl)  minutes 
schedule  a  east. unit. generator  in 

gamma. f(r3"alphal,betal,sl)  minutes 
print  4  lines  thus 

6  event  routine  ■  ■  


Movements    Queue   Next  route 
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Event   Time    U    (N,  S,  E)    (N,S,E)    (N,S,E) 

reserve  arrival,  time.  of.  north,  iinit , 

departure,  time.  of.  north,  unit, 
end.  time.  of.  north,  service, 
waiting,  time.  of.  north,  unit 
and  north. queue. Iength2  as  2000 

reserve  north. queue,  lengthl, 
south. queue,  lengthl, 
east. queue. lengthl  as  2000 

reserve  arrival,  time.  of.  south,  unit, 

departure,  time.  of.  south,  unit, 
end.  time.  of.  south,  service, 
waiting,  time.  of.  south,  unit 
and  south. queue. Iength2  as  2000 

reserve  arrival,  time.  of.  east,  unit, 

departure,  time.  of.  east,  unit, 
end.  time.  of.  east,  service, 
waiting,  time.  of.  east,  unit 
and  east. queue. Iength2  as  2000 

reserve  north. dir, south. dir, east, dir  as  2000 

reserve  AQ.  of.  north,  unit ,AQ.  of.  south.  unit,AQ.  of.  east,  unit, 

AW.  of.  north,  unit, AW.  of.  south,  unit , AW.  of.  east,  unit  as  2000 

END   "INITIALIZE 


EVENT  NORTH.  UNIT. GENERATOR  saving  the  event  notice 
if  jn  <  Nl  or  js  <  Nl  or  je  <  Nl 

schedule  a  new. north. unit  now 
reschedule  this  north. unit,  generator  in 

gamma. f(rl"alphal,  betal,sl)  minutes 
always 
END  '  ■*  NORTH.  UNIT.  GENERATOR 


EVENT  SOUTH.  UNIT.  GENERATOR  saving  the  event  notice 
if  jn  <  Nl  or  js  <  Nl  or  je  <  Nl 

schedule  a  new. south,  unit  now 
reschedule  this  south. unit. generator  in 

gamma.  fCrl'^f alpha  1,  betal,sl)  minutes 
always 
END  ' ' SOUTH.  UNIT.  GENERATOR 


EVENT  EAST. UNIT. GENERATOR  saving  the  event  notice 
if  jn  <  Nl  or  js  <  Nl  or  je  <  Nl 

schedule  a  new. east. unit  now 
reschedule  this  east. unit. generator  in 

gamma.  f(r3'<'alphal,betal,sl)  minutes 
always 
END  '  ■*  EAST.  UNIT.  GENERATOR 
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EVENT  NfEW.  NORTH.  UNIT 
define  U  as  an  real  variable 
create  an  unit 
let   in  =  in  +  1 

current,  time  =  time.  v*hours.  v*niinutes.  v 
arrival,  time.  of.  north.  unit(  in)  =  current,  time 
U  =  uniform. f(0, l,s3) 

if  north,  run  <>  0  or  n. north. queue  >  0 
go  to  'caseS' 
else 
if  U  <  PI 

next. north. route  =  1 
go  to  'easel' 
else 

next. north. route  =  2 
go  to  'case2' 

' easel* 
if  east. run  <>  1 

north. run  =  1 

next. north. route  =  0 

jn  =  jn  +  1 

schedule  an  north. service  giving  unit 

in  gamma.  f(alpha2 ,beta2,s2)  minutes 

departure. time. of . north. unit( jn)  =  current. time 

AQ.  of.  north,  unit(jn)  =  AQ.  north,  unit 
else 
'cases' 

file  this  unit  in  the  north, queue 

north,  queue. Iength2( in)  =  n.  north. queue 
always 
go  to  'end' 

case2 

if  south,  run  =  0  and  east,  run  <>  1 

let  north,  run  =  2 

jn  =  jn  +  1 

next,  north. route  =  0 

schedule  an  north. service  giving  unit 

in  gamma.  f(alpha2,beta2,s2)  minutes 

departure. time. of . north. unit( jn)  =  current. time 

AQ.  of.  north,  unit  (jn)  =  AQ.  north,  vmit 
else 

go  to  'caseS' 
always 

'end' 
print  1  line  with  current. time, U, 

north. run, south. run, east,  run, 
n. north. queue, n. south. queue, n.  east,  queue, 
next. north. route, next. south. route, 
next. east. route  thus 
NNU   **. **   . **       **  **  **       **  **  **       **  **  ** 

END    "NEW.  NORTH.  UNIT 
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EVENT  NEW. SOUTH.  UNIT 
define  U  as  real  variable 
create  an  unit 
is  =  is  +  1 

current. time  =  time. v*hours.  v*minutes.  v 
arrival,  time.  of.  south.  unit(  is)  =  current,  time 
let  U  =  uniform.  f(0, l,s3) 

if  n. south. queue  <>  0  or  south. run  <>  0 
go  to  'cases' 
else 
if  U  <  PI 

next. south. route  =  1 
go  to  'easel' 
else 

next. south. route  =  2 
go  to  'case2' 

easel 

if  north. run  <>  2  and  east,  run  =  0 

south. run  =  1 

next. south,  route  =  0 

js  =  js  +  1 

schedule  an  south. service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 
departure. time. of . south. unit(js)  =  current. time 
AQ.  of.  south,  unit(js)  =  AQ.  south,  unit 
else 

case3 

file  this  unit  in  the  south. queue 
south. queue. Iength2(is)  =  n. south. queue 

always 

go  to  'end' 

case2 

if  north. run  <>  2 

let  south. run  =  2 

js  =  js  +  1 

next. south. route  =  0 

schedule  an  south. service  giving  unit 

in  gamma.  f(alpha2,beta2,s2)  minutes 

departure. time. of. south. unit(js)  =  current. time 

AQ.  of .  south,  unit(js)  =  AQ.  south,  unit 
else 

go  to  'cases' 
always 

'end' 
print  1  line  with  current. time, U, 

north. run, south. run, east,  run, 
n. north. queue, n. south,  queue, n.  east,  queue, 
next. north. route, next,  south,  route, 
next. east,  route  thus 

NSU  VfA'^  VrVf  ^  VfVf  •sV-'c    VcVc    VfV-  VfVf    icic    "k-k  ^Vc    VcVc    -kif 
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END  ' ' NEW.  SOUTH.  UNIT 


EVENT  NEW.  EAST.  UNIT 
define  U  as  real  variable 
create  an  unit 
ie  =  ie  +  1 

current,  time  =  time.  v*hours. v*minutes. v 
arrival,  time.  of.  east,  unit(ie)  =  current,  time 
U  =  uniform.  f(0, 1 ,s3) 

if  n. east. queue  >  0  or  east. rvm  <>  0 
go  to  'caseS' 
else 
if  U  <  PI 

next. east,  route  =  1 
go  to  'easel' 
else 

next. east. route  =  2 
go  to  'case2' 
easel 

if  south. run  <>  1  and  north. run  =  0 
east. run  =  1 
next. east,  route  =  0 
je  =  je  +  1 
schedule  an  east,  service  giving  vinit 

in  gamma. f (alpha2,beta2,s2)  minutes 
departure. time. of . east. unit(je)  =  current. time 
AQ.  of.  east,  unite  je)  =  AQ.  east,  unit 
else 
caseS 

file  this  unit  in  the  east. queue 
east. queue. Iength2( ie)  =  n. east. queue 
always 
go  to  'end' 

'case2' 
if  south,  run  <>  1 
east. run  =  2 
je  =  je  +  1 
next. east. route  =  0 
schedule  an  east. service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 
departure. time. of . east,  unit(je)  =  current. time 
AQ.  of .  east.  unit(  je)  =  AQ.  east,  unit 
else 

go  to  'caseS' 
always 

'end' 
print  1  line  with  current. time, U, 

north. run, south,  run, east,  run, 
n. north. queue, n. south. queue, n. east,  queue, 
next. north. route, next,  south,  route, 
next. east. route  thus 
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NEU   **. **        . **   **  **  **   **  **  **       **  **  ** 


»  I 


END    "NEW.  EAST.  UNIT 


EVENT  NORTH. SERVICE  given  unit 
define  unit  as  a  pointer  variable 
define  U  as  real  variable 
destroy  this  unit 
kn  =  kn  +  1 

current. time  =  time. v*hours. v*minutes. v 
end.  time.  of.  north.  serviceCkn)  =  current,  time 
U  =  uniform. f(0, l,s3) 

if  n. north. queue  >  0  and  U  <  PI 

next. north,  route  =  1 
always 

if  n. north. queue  >  0  and  U  >=  PI 

next. north,  route  =  2 
always 

if  north,  run  =  2 

go  to  'service2' 
else 

north. dir(kn)  =  1 

north. run  =  0 

if  south. run  <>  1  and  next. east,  route  =  1 

east. run  =  1 

je  =  je  +  1 

departure,  time.  of.  east,  unit(je)  =  current,  time 

waiting,  time.  of.  east,  unit(je)  = 

current,  time  -  arrival,  time.  of.  east,  unit(je) 

east. queue. lengthl(je)  =  n. east. queue 

remove  the  first  unit  from  the  east. queue 

AQ.  of.  east.  unit(  je)  =  AQ.  east,  unit 

next. east,  route  =  0 

schedule  an  east. service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 
always 

next 

if  east. run  <>  1  and  next. north. route  =  1 

jn  =  jn  +  1 

north. run  =  1 

departure. time,  of .  north,  unit(jn)  =  current,  time 

waiting,  time,  of .  north,  unit(jn)   = 

current. time  -  arrival. time. of . north. unit(jn) 

north. queue. lengthl(jn)  =  n. north. queue 

remove  the  first  unit  from  the  north. queue 

AQ.  of.  north,  unit(jn)  =  AQ.  north,  unit 

next. north,  route  =  0 

schedule  an  north. service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 
always 
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if  east. run  <>  1  and  next. north. route  =  2  and  south,  run  =  0 

jn  =  jn  +  1 

north. run  =  2 

departure. time. of . north. unit(jn)  =  current. time 

waiting,  time.  of.  north,  unit(jn)   = 

current,  time  -  arrival,  time,  of .  north,  unit(jn) 

north,  queue,  lengthl(jn)   =  n.  north,  queue 

remove  the  first  unit  from  the  north. queue 

AQ.  of.  north,  unit(jn)  =  AQ.  north,  unit 

next. north,  route  =  0 

schedule  an  north. service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 
always 
go  to  'print' 

' service2' 
north. dir(kn)  =  2 
north. run     =  0 

if  next. east,  route  =  1 

je  =  je  +  1 

east. run  =  1 

departure,  time.  of.  east,  unit(je)  =  current,  time 

waiting,  time.  of.  east,  unit(je)  = 

current. time  -  arrival. time. of . east,  unit(je) 

east. queue. lengthl(je)  =  n. east. queue 

remove  the  first  unit  from  the  east. queue 

AQ.  of .  east.  unit(  je)  =  AQ.  east,  unit 

next. east. route  =  0 

schedule  an  east. service  giving  unit 

in  gamma.  f(alpha2,beta2,s2)  minutes 
always 

if  east. run  =  0  and  next. south. route  =  1 

js  =  js  +  1 

south,  run  =  1 

departure. time. of. south,  unit(js)  =  current,  time 

waiting,  time.  of.  south,  unit(js)  =  current,  time  - 

arrival,  time.  of.  south.  unit(  js) 

south. queue. lengthl(js)  =  n. south. queue 

remove  the  first  unit  from  the  south. queue 

AQ.  of.  south,  unit(js)  =  AQ.  south,  unit 

next,  south,  route  =  0 

schedule  an  south. service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 
always 
if  next,  south,  route  =  2 

js  =  js  +  1 

south. run  =  2 

departure. time. of .  south,  unit(js)  =  current,  time 

waiting,  time,  of .  south,  unit(js)  = 

current. time  -  arrival. time. of . south. unit(js) 

south. queue. lengthl(js)  =  n. south. queue 

remove  the  first  unit  from  the  south, queue 

AQ.  of.  south,  unit(js)  =  AQ.  south,  unit 

next. south. route  =  0 
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schedule  an  south. service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 
always 


go  to  'next' 


'print' 


print  1  line  with  current,  time, U, 

north. run, south,  run, east,  run, 
n. north. queue, n.  south,  queue, n.  east,  queue, 
next. north,  route, next,  south,  route, 
next. east,  route  thus 
NS    **.  **   .  **   **  '^''^  ^'^   **  '^*  **   **  **  ** 
END   "NORTH.  SERVICE 


EVENT  SOUTH.  SERVICE  given  unit 
define  U  as  real  variable 
define  unit  as  a  pointer  variable 
ks  =  ks  +  1 
destroy  this  unit 

current. time  =  time. v*hours.  v^minutes.  v 
end.  time.  of.  south,  service(ks)  =  current,  time 
U  =  uniform.  f(0, l,s3) 

if  n. south. queue  >  0  and  U  <  PI 

next. south,  route  =  1 
always 

if  n. south. queue  >  0  and  U  >=  PI 

next. south,  route  =  2 
always 

if  south. run  =  2 

go  to  'service2' 
else 

south,  dir(ks)  =  1 

south. run     =  0 

if  next. east. route  =  1  and  north. run  =  0 

je  =  je  +  1 

east. run  =  1 

departure,  time.  of.  east,  unit(je)  =  current,  time 

waiting,  time,  of .  east,  unit(je)  = 

current,  time  -  arrival,  time,  of .  east,  unit(je) 

east. queue. lengthl( je)  =  n. east. queue 

remove  the  first  unit  from  the  east. queue 

AQ.  of .  east.  unit(  je)  =  AQ.  east,  unit 

next. east,  route  =  0 

schedule  an  east. service  giving  unit 

in  gamma.  f(alpha2,beta2,s2)  minutes 
always 

if  next. east,  route  =  2 
je  =  je  +  1 
east. run  =  2 
departure,  tire.  of.  east,  unit(je)  =  current,  time 
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waiting,  time. of . east. unit(je)  =  current. time  - 

arrival,  time.  of.  east.  unit(  je) 
east. queue. lengthl( je)  =  n. east. queue 
remove  the  first  unit  from  the  east. queue 
AQ.  of.  east.  unit(  je)  =  AQ.  east,  unit 
next. east,  route  =  0 
schedule  an  east. service  giving  unit 

in  gamma.  f(alpha2,beta2,s2)  minutes 
always 
'next' 
if  next,  north,  route  =  2   and  east,  run  <>  1 
jn  =  jn  +  1 
north. run  =  2 

departure,  time. of . north, unit(jn)  =  current. time 
waiting,  time,  of .  north,  unit(jn)  = 

current. time  -  arrival. time. of . north,  unit(jn) 
north,  queue,  lengthl(jn)  =  n.  north,  queue 
remove  the  first  unit  from  the  north. queue 
AQ.  of .  north,  unit(jn)  =  AQ.  north,  unit 
next,  north,  route  =  0 
schedule  an  north. service  giving  unit 

in  gamma.  f(alpha2,beta2,s2)  minutes 
always 

if  east. run  =  0  and  north. run  <>  2  and  next. south. route 

js  =  js  +  1 

south,  run  =  1 

departure. time.  of.  south,  unit(js)  =  current,  time 

waiting,  time.  of.  south,  unit(js)  = 

current. time  -  arrival. time,  of .  south,  unit(js) 

south. queue. lengthl(js)  =  n. south. queue 

remove  the  first  unit  from  the  south. queue 

AQ.  of .  south,  unit(js)  =  AQ.  south,  unit 

next. south. route  =  0 

schedule  an  south,  service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 
always 


=  1 


if  north,  run  <>  2  and  next. south. route  =  2 

js  =  js  +  1 

south. run  =  2 

departure,  time.  of.  south,  unit(js)  =  current,  time 

waiting,  time.  of.  south,  unit(js)  = 

current,  time  -  arrival,  time,  of .  south,  unit(js) 

south. queue. lengthl(js)  =  n. south. queue 

remove  the  first  unit  from  the  south. queue 

AQ.  of.  south,  unit(js)  =  AQ,  south,  unit 

next, south, route  =  0 

schedule  an  south. service  giving  unit 

in  gamma,  f (alpha2,beta2,s2)  minutes 
always 
go  to  'print' 

' service2' 
south. dir(ks)  =  2 
south, run     =  0 
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go  to  'next' 


print 


print  1  line  with  current. time, U, 

north. run, south,  run, east,  run, 
n. north. queue, n. south,  queue, n.  east,  queue, 
next. north. route, next,  south,  route, 
next. east. route  thus 

END  ' ' SOUTH.  SERVICE 


EVENT  EAST. SERVICE  given  unit 
define  U  as  real  variable 
define  unit  as  a  pointer  variable 
destroy  this  unit 
ke  =  ke  +  1 

current. time  =  time. v*hours. v*minutes.  v 
end.  time.  of.  east,  service(ke)  =  current,  time 
U  =  uniform.  f(0, l,s3) 

if  the  east. queue  is  not  empty  and  U  <  PI 
next. east. route  =  1 

always 

if  the  east. queue  is  not  empty  and  U  >=  PI 

next. east,  route  =  2 
always 

if  east. run  =  2 

go  to  'service2' 
else 

east. dir(ke)  =  1 

east, run     =  0 

if  next. north,  route  =  1 

jn  =  jn  +  1 

north. run  =  1 

departure. time. of . north. unit(jn)  =  current,  time 

waiting,  time,  of .  north,  unit(jn)  = 

current. time  -  arrival. time. of . north,  unit(jn) 

north. queue. lengthl(jn)  =  n. north. queue 

remove  the  first  unit  from  the  north. queue 

AQ.  of.  north,  unit(jn)  =  AQ.  north,  unit 

next. north,  route  =  0 

schedule  an  north. service  giving  unit 

in  gamma.  f(alpha2,beta2,s2)  minutes 
always 

if  next,  north. route  =  2  and  south,  run  <>  2 
jn  =  jn  +  1 
north. run  =  2 

departure. time,  of .  north,  unit(jn)  =  current. time 
waiting,  time.  of.  north,  unit(jn)  = 

current,  time  -  arrival,  time.  of.  north,  unit(jn) 
north. queue. lengthl(jn)  =  n. north. queue 
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remove  the  first  unit  from  the  north. queue 
AQ.  of .  north,  unit(jn)  =  AQ.  north,  unit 
next. north,  route  =  0 
schedule  an  north. service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 
always 

'next* 
if  next,  south,  route  =  1  and  north. run  <>  2 
js  =  js  +  1 
south. run  =  1 

departure. time. of. south,  unit(js)  =  current. time 
waiting,  time,  of .  south,  unit(js)  = 

current,  time  -  arrival,  time,  of .  south,  unit(js) 
south. queue. lengthl(js)  =  n. south. queue 
remove  the  first  unit  from  the  south,  queue 
AQ. of . south. unit(js)  =  AQ. south. unit 
next. south,  route  =  0 
schedule  an  south. service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 
always 

if  north,  run  =  0  and  south,  run  <>  1  and  next. east. route  =  1 

je  =  je  +  1 

east. run  =1 

departure,  time.  of.  east,  unit(je)  =  current,  time 

waiting,  time,  of .  east.  unit(  je)  = 

current. time  -  arrival. time. of . east. unit(je) 

east,  queue,  lengthl(je)  =  n.  east,  queue 

remove  the  first  unit  from  the  east. queue 

AQ.  of.  east.  unit(  je)  =  AQ.  east,  unit 

next. east,  route  =  0 

schedule  an  east. service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 
always 

if  south,  run  <>  1  and  next. east. route  =  2 

je  =  je  +  1 

east,  run  =  2 

departure,  time.  of.  east,  unit(je)  =  current,  time 

waiting,  time,  of .  east,  unit(je)  = 

current,  time  -  arrival,  time.  of.  east,  unit(je) 

east. queue,  lengthl(je)  =  n. east. queue 

remove  the  first  unit  from  the  east. queue 

AQ.  of .  east.  unit(  je)  =  AQ.  east,  unit 

next. east,  route  =  0 

schedule  an  east. service  giving  unit 

in  gamma. f (alpha2,beta2,s2)  minutes 
always 
go  to  'print' 

*service2' 
east. dir(ke)  =  2 
east. run  =  0 

go  to  'next' 
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'print' 
print  1  line  with  current. time, U, 

north. run, south. run, east,  run, 
n. north. queue, n. south. queue, n. east,  queue, 
next. north. route, next,  south,  route, 
next. east. route  thus 

£Jg  icif^  **  ^  VfVf  rlfk    id;    iii(  icit    **    **  Vr*    ici(    icic 

END  "EAST.  SERVICE 


ROUTINE  to  SINGLE. RUN.  REPORT 
define  nn,ns,ne  as  integer  variables 
print  3  line  thus 

Single. run. report 


NO    AT     Ql    WT      ST      ET     Q2   DIR 

for  nn  =  1  to  in 
do 

print  1  line  with  nn, 
arrival,  time.  of.  north,  unit(nn) , 
north. queue. Iength2(nn) , 
waiting,  time.  of.  north,  unit(nn) , 
departure,  time.  of.  north,  unit(nn) , 
end.  time.  of.  north,  service(nn) , 
north. queue,  lengthl(nn) , north,  dir(nn)  thus 

*\  *r**  f\  *\      *s  t\  f\  *•  *'.  **      *\  **  **  *\      *««»  «v  *\      •»75^  7r#r  tstt 

loop 

for  ns  =  1  to  is 
do 

print  1  line  with  ns , 
arrival,  time.  of.  south,  unit(ns) , 
south. queue. Iength2(ns) , 
waiting,  time.  of.  south,  unit(ns) , 
departure,  time.  of.  south,  unit(ns) , 
end.  time.  of.  south,  service(ns) , 
south. queue,  lengthl(ns) , south,  dir(ns)  thus 

#*/**»  **   *»  **   *»  7T**  *»   **  *%   r\  •V?*  *»   *S  4\4\  #%  ««  #%  #«  #*7* 

loop 

for  ne  =   1  to   ie 
do 

print  1  line  with  ne, 
arrival,  time.  of.  east,  unit(ne) , 
east. queue. Iength2(ne) , 
waiting,  time,  of .  east,  unit(ne) , 
departure,  time.  of.  east,  unit(ne)  , 
end.  time.  of.  east,  service(ne) , 
east. queue. lengthl(ne) ,east.  dir(ne)  thus 

iticic  VfA^  VrVf  Vr*  *Vf_  i;-h  id;  ^  ifi;  **_  **  **  ** 

loop 

END  ' ' SINGLE.  RUN.  REPORT 
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ROUTINE  to  SUM. REPORT 
define  ii  as  integer  variables 
define  sura. of .  north,  wt, 
sum.  of.  south,  wt, 

sum. of. east, wt  as  1-dimensional  real  array 
define  sum. north,  wt , 
sum.  south,  wt , 

sum. east. wt  as  real  variables 
reserve  sum. of . north. wt, 
sum.  of.  south,  wt , 
sum.  of.  east,  wt  as  2000 

for  ii  =  1  to  Nl 

do 

sum.  north,  wt  =  sura,  north,  wt  + 

waiting,  tirae.  of.  north.  unit(  ii) 
sum.  of.  north,  wt(ii)  =  sum.  north,  wt 
AW.  of.  north,  unit(ii)  =  sura.  of.  north.  wt(ii)/ii 

sum.  south,  wt  =  sum.  south,  wt  + 

waiting,  time.  of.  south.  unit(  ii) 
sum.  of.  south.  wt(  ii)  =  sum.  south,  wt 
AW.  of.  south,  unit(ii)  =  sum.  of.  south.  wt(  ii)/ii 

sum.  east,  wt  =  sum.  east,  wt  + 

waiting,  time.  of.  east.  unit(  ii) 
sum.  of.  east.  wt(  ii)  =  sura.  east,  wt 
AW.  of.  east.  unit(  ii)  =  sura.  of.  east.  wt(ii)/ii 

' 'print  running  average  of  waiting  tirae  and  queue  length, 
if  ii  =  1 
print  3  line  thus 
========   Sura,  report  routine   ====== 

No    AQn    AWn    AQs    AWs    AQe    AWe 


always 

print  1  line  with  ii, 

AQ.  of.  north.  unit(  ii)  ,  AW.  of.  north.  unit(  ii)  , 
AQ.  of.  south.  unit(  ii)  , AW.  of.  south.  unit(  ii)  , 
AQ.  of.  east.  unit(  ii) , 
AW.  of.  east.  unit(  ii)  thus 

loop 

release  sum.  of.  north,  wt  ,sura.  of.  south.  wt,sum.  of.  east,  wt 
stop 
END  ^ ' SUM. REPORT 
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APPENDIX  C.      A  PROGRAM  LISTING  OF  MODEL  2. 

EVENT  NEW. NORTH.  UNIT 
define  U  as  an  real  variable 
create  an  unit 
let   in  =  in  +  1 

current. time  =  time.  V^hours. v*minutes. v 
arrival. time. of . north. unit( in)  =  current. time 
U  =  uniform. f(0 J l,s3) 

if  north. run  <>  0  or  n. north. queue  >  0 
go  to  'caseS' 
else 
if  U  <  PI 

next. north. route  =  1 
go  to  'easel' 
else 

next. north. route  =  2 
go  to  'case2' 

'     1  t 
easel 

if  east,  run  <>  1 

north. run  =  1 

next,  north,  route  =  0 

jn  =  jn  +  1 

schedule  an  north. service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 

departure. time. of . north. unit(jn)  =  current. time 

AQ.  of.  north,  unit(jn)  =  AQ.  north,  unit 
else 
' case3' 

file  this  unit  in  the  north. queue 

north. queue. Iength2( in)  =  n. north. queue 
always 
go  to  'end' 

'case2' 
if  south. run  <>  2 

let  north,  run  =  2 

jn  =  jn  +  1 

next. north. route  =  0 

schedule  an  north. service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 

departure,  time.  of.  north,  unit(jn)  =  current,  time 

AQ.  of.  north,  unit(jn)  =  AQ.  north,  unit 
else 

go  to  'cases' 
always 

'end' 
print  1  line  with  current. time, U, 

north. run, south,  run, east,  run, 

n. north. queue, n. south. queue, n.  east,  queue. 
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next. north. route, next. south. route, 
next. east. route  thus 


END    "NEW.  NORTH.  UNIT 


EVENT  NEW. SOUTH. UNIT 
define  U  as  real  variable 
create  an  unit 
is  =  is  +  1 

current. time  =  time. v*hours. v*minutes.  v 
arrival,  time.  of.  south.  unit(  is)  =  current,  time 
let  U  =  uniform. f( 0,1, s3) 

if  n. south. queue  <>  0  or  south. run  <>  0 
go  to  'caseS' 
else 
if  U  <  PI 

next. south. route  =  1 
go  to  'easel' 
else 

next. south. route  =  2 
go  to  'case2' 

easel 

if  east,  run  <>  2 

south. run  =  1 

next. south. route  =  0 

js  =  js  +  1 

schedule  an  south,  service  giving  unit 

in  gamma. f(alpha2,beta2 ,s2)  minutes 
departure. time. of . south. unit(js)  =  current. time 
AQ.  of .  south,  unit(js)  =  AQ.  south,  unit 
else 

case3 

file  this  unit  in  the  south. queue 
south. queue.  Iength2(is)  =  n. south,  queue 

always 

go  to  'end' 

'case2' 
if  north. run  <>  2 

let  south,  run  =  2 

js  =  js  +  1 

next. south. route  =  0 

schedule  an  south. service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 
departure,  time.  of.  south,  unit(js)  =  current,  time 
AQ.  of.  south,  unit(js)  =  AQ.  south,  unit 


else 

go 
always 


to  'cases' 
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'end' 
print  1  line  with  current. time, U, 

north. run  J  south,  run, east,  run, 
n. north. queue, n. south,  queue, n.  east,  queue, 
next. north. route, next,  south,  route, 
next. east,  route  thus 
NSU  **.  **   .  **   **  **  **   **  **  **   *'^  **  ** 

END  ' 'NEW. SOUTH.  UNIT 


EVENT  NEW. EAST.  UNIT 
define  U  as  real  variable 
create  an  unit 
ie  =  ie  +  1 

current. time  =  time. v^hours. V"minutes. v 
arrival. time. of . east. unit( ie)  =  current. time 
U  =  uniform.  f(0, l,s3) 

if  n. east. queue  >  0  or  east. run  <>  0 
go  to  'caseS' 
else 
if  U  <  PI 

next. east. route  =  1 
go  to  'easel' 
else 

next. east. route  =  2 
go  to  'case2' 
' easel' 
if  north,  run  <>  1 
east. run  =  1 
next. east,  route  =  0 
je  =  je  +  1 
schedule  an  east. service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 
departure. time. of . east. unit( je)  =  current. time 
AQ.  of.  east.  unit(  je)  =  AQ.  east,  unit 
else 
case3 

file  this  unit  in  the  east. queue 
east. queue. Iength2( ie)  =  n. east,  queue 
always 
go  to  'end' 

case2 

if  south. run  <>  1 
east. run  =  2 
je  =  je  +  1 
next. east,  route  =  0 
schedule  an  east. service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 
departure. time. of . east,  unit(je)  =  current,  time 
AQ.  of.  east.  unit(  je)  =  AQ.  east,  unit 
else 

go  to  'caseS' 
always 
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'end' 
print  1  line  with  current. time, U, 

north. run, south. run, east,  run, 
n. north. queue, n. south. queue, n.  east,  queue, 
next. north. route, next,  south,  route, 
next. east,  route  thus 
NEU  **.  **   .  **   **  **  **   '''■*  ''^'^  **   **  **  ** 
END  "NEW.  EAST.  UNIT 


EVENT  NORTH.  SERVICE  given  unit 
define  unit  as  a  pointer  variable 
define  U  as  real  variable 
destroy  this  unit 
kn  =  kn  +  1 

current. time  =  time.  v*hours. v*minutes. v 
end.  time.  of.  north,  service(kn)  =  current,  time 
U  =  uniform. f(0, l,s3) 

if  n. north,  queue  >  0  and  U  <  PI 

next. north. route  =  1 
always 

if  n. north,  queue  >  0  and  U  >=  PI 

next. north. route  =  2 
always 

if  north,  run  =  2 

go  to  'service2' 
else 

north. dir(kn)  =  1 

north. run  =  0 

if  next. east. route  =  1 

east. run  =  1 

je  =  je  +  1 

departure. time. of . east. unit(je)  =  current. time 

waiting,  time.  of.  east,  unit(je)  = 

current. time  -  arrival. time. of . east,  unit(je) 

east. queue,  lengthl(je)  =  n. east,  queue 

remove  the  first  unit  from  the  east. queue 

AQ.  of .  east.  unit(  je)  =  AQ.  east,  unit 

next. east,  route  =  0 

schedule  an  east. service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 
always 

'next* 
if  east. run  <>  1  and  next. north. route  =  1  . 
jn  =  jn  +  1 
north. run  =  1 

departure,  time. of . north. unit(jn)  =  current. time 
waiting,  time,  of .  north,  unit(jn)   = 

current,  time  -  arrival,  time.  of.  north,  unit(jn) 
north. queue,  lengthl(jn)  =  n. north,  queue 
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remove  the  first  unit  from  the  north. queue 
AQ.  of.  north,  unit(jn)  =  AQ.  north,  unit 
next. north,  route  =  0 
schedule  an  north. service  giving  unit 

in  gamma. f (alpha2 jbeta2,s2)  minutes 
always 

if  next. north. route  =  2  and  south,  run  <>  2 

jn  =  jn  +  1 

north,  run  =  2 

departure,  time,  of .  north,  unit(jn)  =  current,  time 

waiting,  time.  of.  north,  unit(jn)  = 

current. time  -  arrival. time,  of .  north,  unit(jn) 

north. queue,  lengthl(jn)  =  n. north. queue 

remove  the  first  unit  from  the  north,  queue 

AQ.  of.  north,  unit(jn)  =  AQ.  north,  unit 

next. north,  route  =  0 

schedule  an  north. service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 
always 
go  to  'print' 

service2 

north. dir(kn)  =  2 

north,  run     =  0 

if  next. south,  route  =  2 

js  =  js  +  1 

south. run  =  2 

departure,  time.  of.  south,  unit(js)  =  current,  time 

waiting,  time.  of.  south,  unit(js)  = 

current. time  -  arrival. time. of . south. unit(js) 

south. queue. lengthl( js)  =  n. south. queue 

remove  the  first  unit  from  the  south. queue 

AQ.  of .  south.  unit(  js)  =  AQ.  south,  unit 

next. south,  route  =  0 

schedule  an  south. service  giving  unit 

in  gamma.  f(alpha2,beta2,s2)  minutes 
always 
go  to  'next' 

print 

print  1  line  with  current. time, U, 

north,  run, south,  run, east. run, 
n.  north,  queue, n.  south,  queue, n.  east,  queue, 
next. north,  route, next,  south,  route, 
next. east. route  thus 
fi^S   **.  **    .  **   ''f*  '''^'^   '^*   **  *'"''  *''<■   **  **  ^'^ 
END   "NORTH.  SERVICE 


EVENT  SOUTH. SERVICE  given  unit 
define  U  as  real  variable 
define  unit  as  a  pointer  variable 
ks  =  ks  +  1 
destroy  this  unit 
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current. time  =  time. v*hours. v*minutes.  v 
end.  time.  of.  south,  service(ks)  =  current,  time 
U  =  uniform.  f( 0, 1 ,s3) 

if  n. south. queue  >  0  and  U  <  PI 

next. south,  route  =  1 
always 

if  n. south. queue  >  0  and  U  >=  PI 

next. south,  route  =  2 
always 

if  south,  run  =  2 

go  to  'service2' 
else 

south. dir(ks)  =  1 

south,  run     =  0 

if  next. east,  route  =  2 

je  =  je  +  1 

east. run  =  2 

departure. time.  of.  east,  unit(je)  =  current,  time 

waiting,  time.  of.  east,  unit(je)  =  current,  time  - 

arrival,  time.  of.  east.  unit(  je) 

east. queue. lengthl(je)  =  n. east. queue 

remove  the  first  unit  from  the  east. queue 

AQ.  of .  east.  unit(  je)  =  AQ.  east,  unit 

next. east. route  =  0 

schedule  an  east. service  giving  unit 

in  gamma. f (alpha2,beta2,s2)  minutes 
always 

' next ' 
if  east. run  <>  2  and  next,  south,  route  =  1 
js  =  js  +  1 
south. run  =  1 

departure. time. of.  south,  unite js)  =  current,  time 
waiting,  time.  of.  south,  unit(js)  = 

current. time  -  arrival. time. of . south. unit(js) 
south. queue. lengthl(js)  =  n. south. queue 
remove  the  first  unit  from  the  south. queue 
AQ.  of.  south,  unit(js)  =  AQ.  south,  unit 
next. south,  route  =  0 
schedule  an  south. service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 
always 

if  north. run  <>  2  and  next. south. route  =  2 
js  =  js  +  1 
south. run  =  2 

departure. time. of. south,  unit(js)  =  current,  time 
waiting,  time.  of.  south,  unit(js)  = 

current,  time  -  arrival,  time.  of.  south,  vinit(js) 
south. queue. lengthl( js)  =  n. south. queue 
remove  the  first  unit  from  the  south. queue 
AQ.  of.  south,  unit(js)  =  AQ.  south,  unit 
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next. south. route  =  0 

schedule  an  south,  service  giving  unit 

in  gamma. f (alpha2,beta2,s2)  minutes 
always 
go  to  'print' 

serviceZ 

south. dir(ks)  =  2 

south,  run     =  0 


if  next. north. route  =  2  and  south,  run  <>  2 

jn  =  jn  +  1 

north. run  =  2 

departure. time. of.  north,  unit(jn)  =  current,  time 

waiting,  time.  of.  north,  unit(jn)  = 

current,  time  -  arrival,  time.  of.  north,  unit(jn) 

north. queue. lengthl(jn)  =  n. north. queue 

remove  the  first  unit  from  the  north. queue 

AQ.  of.  north,  imit(  jn)  =  AQ.  north,  unit 

next. north,  route  =  0 

schedule  an  north. service  giving  unit 

in  gamma,  f (alpha2 jbeta2,s2)  minutes 
always 
go  to  'next' 

'print' 
print  1  line  with  current,  time, U, 

north. run, south,  run, east,  run, 
n. north. queue, n. south. queue, n.  east,  queue, 
next. north. route, next,  south,  route, 
next, east. route  thus 


SS 


^•jV    Vf* 


*Vf 


icit 


i(it    iti;    ick  i^-h    ich    ick 


END 


SOUTH.  SERVICE 


EVENT  EAST. SERVICE  given  unit 
define  U  as  real  variable 
define  unit  as  a  pointer  variable 
destroy  this  unit 
ke  =  ke  +  1 

current,  time  =  time,  v'^'hours.  v*minutes.  v 
end. time. of . east. service(ke)  =  current. time 
U  =  uniform.  f(0,l,s3) 

if  the  east. queue  is  not  empty  and  U  <  PI 

next. east. route  =  1 
always 

if  the  east. queue  is  not  empty  and  U  >=  PI 

next. east,  route  =  2 
always 

if  east. run  =  2 

go  to  'service2' 
else 

east. dir(ke)  =  1 
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east. run     =  0 

if  next,  north,  route  =  1 

jn  =  jn  +  1 

north. run  =  1 

departure,  time.  of.  north. unit(jn)  =  current. time 

waiting,  time.  of.  north.  unit(  jn)  = 

current. time  -  arrival. time. of . north. unit(jn) 

north,  queue. lengthl( jn)  =  n.  north. queue 

remove  the  first  unit  from  the  north. queue 

AQ.  of.  north.  unit(  jn)  =  AQ.  north,  unit 

next. north,  route  =  0 

schedule  an  north,  service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 
always 

'next' 
if  north,  run  <>  1  and  next. east. route  =  1 
je  =  je  +  1 
east. run  =1 

departure,  time.  of.  east,  unit(je)  =  current,  time 
waiting,  time.  of.  east,  unit(je)  = 

current. time  -  arrival. time. of .  east,  unit(je) 
east. queue. lengthl(je)  =  n. east. queue 
remove  the  first  unit  from  the  east. queue 
AQ.  of.  east.  unit(  je)  =  AQ.  east,  unit 
next. east,  route  =  0 
schedule  an  east. service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 
always 

if  south. run  <>  1  and  next. east. route  =  2 

je  =  je  +  1 

east. run  =  2 

departure. time,  of .  east,  unit(je)  =  current. time 

waiting,  time.  of.  east,  unit(je)  = 

current. time  -  arrival. time,  of . east,  unit(je) 

east,  queue.  lengthl( je)  =  n. east. queue 

remove  the  first  unit  from  the  east. queue 

AQ.  of.  east.  unit(  je)  =  AQ.  east,  unit 

next. east,  route  =  0 

schedule  an  east,  service  giving  unit 

in  gamma.  f(alpha2,beta2,s2)  minutes 
always 
go  to  'print' 

' service2' 
east. dir(ke)  =  2 
east,  run  =  0 

if  next,  south,  route  =  1 
js  =  js  +  1 
south. run  =  1 

departure,  time.  of.  south,  unit(js)  =  current,  time 
waiting,  time.  of.  south,  unit(js)  = 

current,  time  -  arrival,  time.  of.  south,  unit(js) 
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south. queue. lengthl( js)  =  n. south. queue 

remove  the  first  unit  from  the  south. queue 

AQ.  of.  south,  unit(js)  =  AQ.  south,  unit 

next. south. route  =  0 

schedule  an  south. service  giving  unit 

in  gamma. f(alpha2,beta2,s2)  minutes 
always 
go  to  'next' 

'print' 
print  1  line  with  current,  time, U, 

north. run, south. run, east. run, 
n.  north,  queue, n.  south,  queue, n.  east,  queue, 
next. north. route, next. south,  route, 
next. east. route  thus 
ES   **.  **   .  **   *''^  *'^  **   **  '^''^   **   **  '^'^   ** 

END  ' 'EAST.  SERVICE 
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